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(57) Abstract 

A system and method for estimating Doppler shift compensation in a mobile communication system is described. A channel impulse 
response determination circuit determines a channel impulse response for the channel on which the signal is received. A reference vector is 
generated using a set of known symbols and the channel impulse response, and an error vector is then generated by comparing the reference 
vector with received samples of the signal. In addition, the noise variance of the incoming signal can be determined. Then, a Doppler 
correction factor is estimated having a real part representing amplitude correction and an imaginary part representing phase correction. 
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(57) Abstract 

A system and method for estimating Doppler shift compensation in a mobile communication system is described. A channel impulse 
response determination circuit determines a channel impulse response for the channel on which the signal is received. A reference vector is 
generated using a set of known symbols and the channel impulse response, and an error vector is then generated by comparing the reference 
vector with received samples of the signal. In addition, the noise variance of the incoming signal can be determined. Then, a Doppler 
correction factor is estimated having a real part representing amplitude correction and an imaginary part representing phase correction. 
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ESTIMATION OF DOPPLER SHIFT COMPENSATION 



IN A MOBILE COMMUNICATION SYSTEM 



The present invention relates to estimating a Doppler shift 
compensation in a mobile communication system. 

In a mobile communication system, signals which are transmitted 
from mobile stations moving relative to a base station are 
subject to the well-known Doppler effect, which causes a 
frequency shift in the frequency received at the base station 
relative to that which was transmitted by the mobile station. 
This frequency shift is referred to herein as the Doppler shift. 
The Doppler shift is dependent upon the speed and direction of 
movement of the mobile station relative to the base station. 
Thus, the Doppler effect can provide an increase or a decrease 
in the frequency, depending on the direction of movement of the 
mobile station relative to the base station. The magnitude of 
the Doppler shift is dependent on the speed with which the mobile 
station is moving relative to the base station. 

Another area of frequency offsets arises when the transmitter 
unit and the receiver unit are wrongly synchronised. 

Existing mobile communication installations provide a form of 
Doppler compensation, in that the frequency detection circuitry 
within the base station which selects a particular signal on a 
particular channel can take into account a certain amount of 
Doppler shift in the signal. 

AU 664626 relates to a method and circuit arrangement for 
compensating for the Doppler shift in a radio signal propagating 
between a base station and a mobile station when the mobile 
station approaches and moves past the base station. As a mobile 
station approaches the base station with a decreasing propagating 
time delay, the propagating time delay is integrated at intervals 
to determine how the propagating time delay varies with time. 
This information is used to effect a change in the radio signal 
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frequency at a particular time to compensate for the sudden 
Doppler shift as the mobile station moves past the base station. 
In a TDMA system, the interrogation occurs at intervals equal to 
an integral number of time frames. In the GSM standard a time 
frame comprises eight consecutive time slots and a single 
transmission burst passes between a particular mobile station and 
base station in any one time frame. 

Thus, in this method of Doppler compensation decisions are made 
in response to past and incoming signals to improve the reception 
at a future time, i.e. it is a reactive system. It would be 
desirable to implement a system which can actively compensate for 
Doppler shifts in an incoming signal in real time. 

Moreover, due to fading a signal subject to the Doppler effect 
may have a variation in amplitude as well as frequency. It would 
be desirable to estimate the Doppler effect on amplitude as well 
as frequency. 

According to one aspect of the invention there is provided a 
system for generating a Doppler correction factor for 
compensating for Doppler shift in a signal transmitted between 
a mobile station and a base station in a mobile communication 
system, the system comprising: a channel impulse response 
determination circuit for determining a channel impulse response 
for the channel on which the signal is received; reference 
circuitry for generating a reference vector from a set of known 
symbols and the channel impulse response; means for generating 
an error vector from the reference vector and received samples 
of the signal; means for determining noise variance of the 
incoming signal; and an estimator for generating a Doppler 
correction factor having a real part representing amplitude 
correction and an imaginary part representing phase correction, 
the estimator utilising the noise variance, the error vector, the 
samples of the received signal and the sampling times. 



According to the GSM standard, a TDMA transmission burst 
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comprises a training sequence which is normally used to determine 
the channel impulse response for the channel on which the signal 
is received. This is done by a convolution of the received 
training sequence with a stored version of the training sequence. 
Thus, the training sequence can constitute the set of known 
symbols for use in generating the reference vector, 
alternatively, the set of known symbols can comprise estimated 
symbols from the received signal samples. 

In a TDMA system, the signal comprises a sequence of transmission 
bursts in respective time slots. The Doppler correction factor 
can be generated for each transmission burst and used to correct 
signal samples in that burst. 

In addition for correcting for the Doppler effect, the Doppler 
correction factor also corrects for frequency offsets arising 
from wrong frequency synchronisation between transmitter and 
receiver units . 

According to a further aspect of the present invention there is 
provided a method for compensating said Doppler shift in a signal 
transmitted between a mobile station and a base station in a 
mobile communication system, the method comprising: determining 
a channel impulse response for the channel on which the signal 
is received; generating a reference vector from a set of known 
symbols in the channel impulse response; generating an error 
vector from the reference vector and received samples of the 
signals; determining noise variance of the incoming signals; 
estimating a Doppler correction factor having a real part 
representing amplitude correction and an imaginary part 
representing phase correction, the estimating step using the 
noise variance, the error vector, the samples of the received 
signals and the sampling times; and applying the Doppler 
correction factor to the received signal samples thereby to 
compensate for Doppler shift. 

For a better understanding of the present invention and to show 
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how the same may be carried into effect, reference will now be 
made by way of example to the accompanying drawings in which: 

Figure 1 is a diagram of a signal burst in a mobile 
communicat ion system ; 

Figure 2 is a block diagram of circuitry for implementing 
Doppler shift compensation; and 

Figure 3 is a diagram illustrating the Estim (3 block. 

Figure 1 illustrates a normal burst in a mobile communication 
system according to the GSM standard. This figure represents a 
burst received at a base station. For a TDMA system according 
to the GSM standard, mobile stations transmit bursts as modulated 
signals on frequency channels allocated by a base station 
controller. One frequency channel may support up to eight 
bursts, each burst associated with a respective call, where each 
call is allocated a time slot in which to send the burst. 
Further details of a TDMA system according to the GSM standard 
are not described herein because they are known to a person 
skilled in the art. 

The normal burst contains two packets of 58 bits (DATA) 
surrounding a training sequence (TRS) of 26 bits. Three tail 
bits (TS) are added at each end of the normal burst. The 
training sequence (TRS) is a predetermined sequence of bits which 
is sent by the mobile station (MS) and is known at the base 
station controller (BSC) . In Figure 1, it is assumed that a. 
received signal sample y k represents the k th bit of the burst, in 
this instance the start of the training sequence. The training 
sequence is utilised at the base station controller to estimate 
the impulse response of the channel over which the burst is sent, 
and, at least in one embodiment of the invention to estimate a 
Doppler compensation factor. The actual information which is 
transmitted is located in the data bits (DATA) of the burst. 

As explained earlier, the environment through which a signal 
passes from a mobile station to a base station can vary 
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considerably, depending, amongst other things, on the distance 
between the mobile station and the base station, and interference 
caused by buildings and other structures in the area. As a 
result, the signal strength and signal quality of the signal 
received at the base station varies widely. Moreover, for moving 
mobile stations, the signal received by the base station is 
subject to a Doppler shift which should be corrected. 

Figure 2 illustrates a circuit 1 suitable for implementing a 
Doppler compensation in a GSM system. It should be understood 
that the various blocks in Figure 2, although illustrated as 
separate interconnected entities, do not necessarily represent 
separate physical entities, but are intended to represent 
diagrammatically the various steps which are carried out. The 
blocks could be implemented as circuits or a suitably programmed 
microprocessor may effect each of the functions which is 
individually assigned to the blocks. 

An antenna 20 receives signals 11 from the mobile stations MS. 
The antenna 20 is connected to RF circuitry 22. This circuitry 
22 operates on the received burst to downshift the frequency to 
the baseband frequency and to sample the burst to provide from 
the analogue signal digital sampled values y ± where i is the 
sample index. The RF circuitry 22 outputs the sampled values 
Yo- • - y n for each burst together with the sampling times t 0 . . .t n/ 
typically at the expected bit rate of the transmitted signal. 
The output of RF circuitry 22 is supplied along line 24 to a 
channel impulse response (C.I.R.) block 10, to a variance and 
error calculator 16 to filtering and equalisation circuitry 12, 
to an Estim (3 block 3 6 for estimating a Doppler correction factor 
and to transforming circuitry 4 0 to enable the application of the 
Doppler correction factor to the burst. 

A memory 32 holds the training sequence x which is the 
predetermined sequence of bits which is sent by the mobile 
station MS as a training sequence and received at the base 
station and converted into samples y^g = Y k -»y k+25 . The reference 
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training sequence x is supplied to a reference generator 14 and 
to the channel impulse response (C.I.R.) block 10. The reference 
generator 14 also receives the estimated channel impulse response 
h from the channel impulse response block 10. 

The C.I.R. block 10 receives the burst, including the received 
training sequence y TRS and calculates an estimated channel impulse 
response h for example by calculating the cross correlation 
between the received training, sequence y TRS and the known training 
sequence x according to equation 1. There are other types of 
estimation methods which could be used. 

h = xcorr (x,y TRS ) (equation 1) 

It will be appreciated that, prior to effecting the cross 
correlation, the known training sequence x, which is stored in 
digital form, is i,q modulated similarly to the manner in which 
the training sequence has been modulated at the MS for 
transmission, according to the GSM standard. The cross 

correlation is done in a known manner to produce a channel 
impulse response in the form of five tap values (h ( j ) j=0to4 ) . 

As is known, the estimated impulse response h is used to 
calculate the expected estimate of the data in the received 
burst, as though the data has been subject to the same average 
noise . 

The C.I.R. block also generates timing advance information r, 
which is used to determine where in the allocated time slot the 
received burst is located. 

For each burst, the estimated channel impulse response h for that 
burst is calculated by the CIR block 10 and is supplied to 
filtering/equalising circuitry which allows the data, DATA ( r ) , 
in that burst to be recovered. As is known, the 

filtering/equalising circuit 12 receives the channel impulse 
response h and timing information r for the received burst to 
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allow the signal to be demodulated, filtered and decoded, to 
recover the data in a known manner. 

The reference generator 14 produces a reference vector, z, which 
is calculated using the convolution of the impulse response and 
the known training sequence. The reference generator 14 performs 
the following calculation: 

2 = h * x (equation 2) 

where x represents the known symbol values depending o the used 
modulation. 

z k represents the kth sample of the vector z and z denotes a 
diagonal matrix of the vector samples z k , that is, z=xh, where x 
is the diagonal matrix of the known training sequence symbols . 

The reference vector z is supplied from the reference generator 
to the variance and error calculator 16. As described above, the 
variance calculator also receives the samples y it including the 
received training sequence y TOS . The variance calculator 
calculates a variance var (a 2 ) according to the following 
equation : 

25 

var = E(|y ± - z ± | A 2) ) (equation 4) 

i=4 



ref_length 

The term ref_length is a constant representing the length of the 
reference vector, z. This is calculated by multiplying the 
number of samples (22) by the bit separation. 

In equation 4, the values of y A are the sampled values of the 
received training sequence for the burst (y k -»y k+25 ) . 

The variance and error calculator also generates an error vector 
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z-y (i.e. error values z—yi for each sample). 

The error vector z-y and the variance a 2 are supplied to the 
Estim (3 block 36. The Estim jS block 36 also receives the signal 
samples and times y^ti from the RF circuitry 22. The Estim (3 
block can be implemented as a processor executing a program or 
by other suitable means. It operates to generate a Doppler 
compensation factor j8 which is a complex number having real and 
imaginary parts. Thus, 0 may be represented as: 

Qf+ j0t 

where or represents amplitude variations due to fading and <f> is 
the Doppler angle. The Estim jS block 36 also includes a memory 
for holding values required to generate the Doppler correction 
factor p, including E{(3) and E(z-y). These terms are 

explained later with reference to Figure 3 . 

The Doppler correction factor 13 together with the signal samples 
y ± and times ti for each burst are supplied to the transforming 
circuitry 4 0 to implement Doppler correction according to the 
following equation: 

Yn-compensated = Yn (1 " ^^n^ " 

Another implementation of the transforming circuitry 4 0 is to 
separate the amplitude and phase components into separate 
correction steps. The phase can be corrected by multiplying the 
samples with e " 3 * ina9Wt n and the amplitude can be corrected in a 
separate step: 

v = v *e-^ ima 9 Wt 

Jr n-phase -compensated I n w n 

v - v .. (1-real {B) t n ) (equation 5) 

Yn-compensated — Y -phase- compensated v - 1 - \ H / -n' ^ 

According to the above, imag(j3)=<£, and real Thus, the 

received signal is either amplified or alternated during the 
burst depending whether real {(3) is negative or positive. 



WO 00/54431 



PCT/EP99/01618 



9 

This generates a set of compensated signal samples y comp which can 
be fed to the filtering/equalising circuitry 12 to generate 
refined values for the data of the burst. 

As an alternative to the above, in which the training sequence 
x is used to supply known symbols as a reference, bits estimated 
in an equalisation circuit 30 can be used. The equalisation 
circuit 30, for example of Viterbi equaliser, receives the 
filtered, demodulated and equalised signal DATA (r) from the 
filtering and equalisation circuitry 12 . The equalisation 
circuit 3 0 operates on a part of the data sequence data of the 
burst (that part having been derived from Est im. block in Figure 
1) to estimate output bits which were sent from the mobile 
station. This output is referred to herein as estim_bits. The 
equalisation circuit 30 operates to make decisions of the bits 
as in known mobile communication systems and thus it will not be 
described further herein. These estimated bits can be used in 
place of the training sequence to supply known values for 
calculating the reference vector z. The estim_bits are supplied 
to a reference circuit 34 with the channel impulse response h 
which generates an error vector z'-y, where z ' =estim_bits*h . 



The theory underlying the generation of the Doppler correction 
factor (3 will now be outlined. 



The problem of (pure) Doppler frequency estimation can be written 
as : 



Xh=Yr+w, where r is a rotation vector, 
Y is a diagonal sample matrix, and 
w is noise. 





r 

Yo 
0 


0 

Yi 


0 
0 




e J * 0 


Y= 


0 




• • Yn-l 

J 


, r= 


e n-l 
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Xh corresponds to convolution between known symbols (diagonal 
matrix x of training sequence symbols) and impulse response h. 

So, in the Doppler frequency estimation, the Doppler angle <f> is 
to be estimated. This problem is nonlinear and therefore a bit 
tricky to be solved, particularly for implementation in real 
time. A linearisation step can be performed assuming e j, * t 0 «l+ j<£t 0 . 
Due to fading the required Doppler frequency correction cannot 
be assumed to be constant in amplitude. 

If the term j</>t 0 has also a real component, the amplitude 
variations can roughly estimated. Thus, the term j0t o is 
replaced by a more general assumption l+/?t 0 , where (3 is complex 
(or+j^t) . 



The linear problem can be written as (for simplification 
z=Xh=[z 0 z 1 . . .z n -i] T ) • 

Zi-Yi (l+/?ti) +w ±/ 



where p= [yo^y^ . . . y^t^] T and /3 is to be estimated. 

A solution to equation 6 can be found by using an LMMSE estimator 

0=E (/?) +C<, u . y) C Ia _ y) u- y> - 1 U-y-E (z-y) ) 

where E(/3) is an expected value for {3, C fiiz _ y} , C (z . y)u . y) are 
covariance matrixes between the two elements, viz E(/?(z-y) H ) and 
E ( (z-y) (z-y) H and E(z-y) is the expected error in the received 
burst . 



which can be reformatted in matrix notation as 



z-y=p/?+w 



(equation 6) 



In this phase it can be assumed that the noise is white gaussian 
noise having a variance of 6 2 and after some manipulation 
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iS=E(i8)+p H (z-y-E(2-y) )/(6 2 /C^+p H p) , where 
C,=E( (P-E(fi) ) (jS-E(jB) ) H ) 



(equation 7) 
(equation 8) 



p H is the complex conjugate transpose matrix of p. 



Equation 7 is used in the Estim f3 block 3 6 to generate /?. The 
Estim (3 block is illustrated schematically in Figure 3. In 
practice, it can be implemented as software executed by a 
processor with a memory for holding the required data. An 
Estimate_j8 block 50 operates on the data (Y i# ti) received in each 
time slots to generate a value for (3 which is used by the 
transforming circuitry 40. An Acc (3 block 52 accumulates and 
averages (3 values over preceding time slots to generate E(jS). 
The Acc B function can be implemented by an adaptive filter. An 
Acc (z-y) block 54 generates an expected error E(z-y) . To do 
this, it can use values for z-y from preceding time slots, or 
E(/3) from the Acc (3 block 52 to estimate a frequency error. A 
Gen C 0 block 56 produces a value for C^. A number of 
possibilities exist for the parameters of equation 7. 

If E(/?) is set to zero, and 5 2 /C p is also set to zero, equation 
7 resolves to a normal linear regression. 

As mentioned above, it is intended that (3 will be generated for 
each time slot. The expected value E(/3) can be generated by 
calculating the average of estimated (3 values from previous time 
slots. Thus, these values can be retained as they are calculated 
for each time slot, and used to produce an average value held in 
the E(jS) block in the Estim 0 block 36. 

Another simplification is to assume that the real part of E((3) 
is zero. This is a reasonable assumption because, in typical 
fading channel situations, no long time amplitude is expected. 
Thus, it is adequate to perform long time averages only in 
respect of error components of frequency. Also, the real part 
could be assumed to be zero if only phase variations are to be 
corrected. 



WO 00/54431 



PCT/EP99/01618 



12 

It is possible to avoid having to 
E(z-y) by derotating the received 
this purpose, /? can be an earlier 
previous time slot. In this case, 

0=E (fi) +p H (z-y) / OVC^+pap) 



calculate the expected error 
samples yi by e " j * ilBagWt . For 
estimated value of /3 from a 
equation 7 changes to 

(equation 9) 



Thus, this requires a preprocessing step (see 50a in Figure 3) 
in which the sample values y ± . are derotated by e -j* im ag<0> ^ p r i or to 
forming the vector p . 

The value for C(/3) can be generated by the GEN block using 
different approaches. It can be a constant value. I can be 
estimated using the information of previous time slots, using for 
example equation 8. It can be made a function of the expected 
value of j3 and/or the noise variance 6 2 . It can be estimated 
using an adaptive filter, or any combination of the preceding. 

The technique is particularly suitable for short delay spread 
channels and pure frequency errors. 
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CLAIMS: 

1. A system for generating a Doppler correction factor for 
compensating for Doppler shift in a signal transmitted between 
a mobile station and a base station in a mobile communication 
system, the system comprising: 

a channel impulse response determination circuit for 
determining a channel impulse response for the channel on which 
the signal is received; 

reference circuitry for generating a reference vector from 
a set of known symbols and the channel impulse response ; 

means for generating an error vector from the reference 
vector and received samples of the signal; 

means for determining noise variance of the incoming signal; 

and 

an estimator for generating a Doppler correction factor 
having a real part representing amplitude correction and an 
imaginary part representing phase correction, the estimator 
utilising the noise variance, the error vector, the samples of 
the received signal and the sampling times. 

2. A system according to claim 1, wherein the Doppler 
correction factor 0 is generated according to the following 
equation : 

/S-E (jS) +p H (z-y-E (z-y) ) / <6 a /C,+p«p> 
wherein E(/3) is an expected value of (3, 

p H is a complex conjugate transpose matrix of p, where p is the 
vector (y ±/ t i ) T / 

z-y is the error vector, 

E(z-y) is the expected error in the received samples, 
5 2 is the noise variance, and 



WO 00/54431 



PCT/EP99/01618 



14 

is a constant . 

3. A system according to claim 2, wherein E(/3) and 6 2 /Cp are set 
to zero. 

4. A system according to claim 1, 2 or 3, which comprises means 
for applying a Doppler correction to the received signal samples 
using the Doppler correction factor /?. 

5. A system according to any preceding claim wherein the signal 
comprises a transmission burst in a time slot in a TDMA mobile 
communication system, and wherein a new Doppler correction factor 
0 is generated for each time slot. 

6 . A system according to claim 5 , which comprises means for 
estimating E(/8) as an average of generated values for (3 from 
previous time slots. 

7. A system according to claim 5, wherein the real part of E(£) 
*is set to be zero. 

8. A system according to any preceding claim, which comprises 
means for derotating the received signal samples by e " 3 * lina9<E( * ,)t , 
where E((3) is an expected value for /S the derotated received 
signal samples then being used to generate the Doppler correction 
factor £ in the present time slot according to the following 
equation: 

/3=E ( /?) +p H ( z-y) / ( SVC^+p"?) . 

9. A method for compensating said Doppler shift in a signal 
transmitted between a mobile station and a base station in a 
mobile communication system, the method comprising: 

determining a channel impulse response for the channel on 
which the signal is received; 

generating a reference vector from a set of known symbols 
in the channel impulse response; 
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generating an error vector from the reference vector and 
received samples of the signals; 

determining noise variance of the incoming signals; 

estimating a Doppler correction factor having a real part 
representing amplitude correction and an imaginary part 
representing phase correction, the estimating step using the 
noise variance, the error vector, the samples of the received 
signals and the sampling times; and 

applying the Doppler correction factor to the received 
signal samples thereby to compensate for Doppler shift. 

10. A method according to claim 9, wherein the signal comprises 
a sequence of transmission bursts in a TDMA mobile communication 
system, and the Doppler correction factor is generated for each 
transmission burst. 

11. A method according to claim 9 or 10, wherein the Doppler 
correction factor (3 is generated according to the following 
equation: 

jS=E (13) +p H (z-y-E (z-y) ) / (fiV^+p^) 
wherein E(j3) is an expected value of /?, 

p H is a complex conjugate transpose matrix of p, where p is the 
vector (yi/ti) 7 , 

z-y is the error vector, 

E(z-y) is the expected error in the received samples, 
6 2 is the noise variance, and 



Cp is a constant. 
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(Any replacement sheet containing such amendments must be referred to under item 1 and annexed to this 
report.) 
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Claims 




Inventive step (IS) 
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Claims 
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Claims 




Industrial applicability (IA) 
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Claims 
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Reasoned statement under Article 35(2) with regard to novelty, inventive step or industrial applicability; 
citations and explanations supporting such statement 

1 . Reference is made to the following document: 

D1 : WO 98 34357 A (NOKIA TELECOMMUNICATIONS OY ;PIIRAINEN OLLI 
(Fl)) 6 August 1998 (1998-08-06) 

2. The document D1 is regarded as being the closest prior art to the subject-matter 
of claim 1 , and shows (the references in parentheses applying to this document; 
see especially Figure 2 and page 7, line 10 to page 13, line 13): 

A system for generating a Doppler correction factor for compensating for Doppler 
shift in a signal transmitted between a mobile station and a base station in a 
mobile communication system, the system comprising: 

a channel impulse response determination circuit (10) for determining a channel 
impulse response for the channel on which the signal is received; 
reference circuitry (14, 34) for generating a reference vector (ref) from a set of 
known signals and the channel impulse response; 

means (36) for generating an error vector (Ph-diff) from the reference vector (ref) 
and received samples of the signal (r); 

means (16) for determining noise variance (a 2 ) of the incoming signal; 
an estimator (18, 38 or 33) for generating a Doppler correction factor, the 
estimator utilising the noise variance, the error vector (Ph-diff), the samples of the 
received signal (r) and the sampling time (h_time). 

The features of claims 1 and 9 which are not known from document D1 relate to 
the Doppler correction factor having a real part representing amplitude correction 
and an imaginary part representing phase correction. In document D1 (see e.g. 
equations 6, 7 and 9), only the phase change resulting from the Doppler effect is 
calculated and corrected. 

Therefore, the subject-matter of claims 1 and 9 is new (Article 33(2) PCT). 
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3. It also involves an inventive step as required by Article 33(3) PCT, because none 
of the available prior art documents suggests to take into account the amplitude 
variation in the signal due to the Doppler effect. 



Re Item VII 

Certain defects in the international application 

1 . The features of the claims are not provided with reference signs placed in 
parentheses (Rule 6.2(b) PCT). 

2. Contrary to the requirements of Rule 5.1 (a)(ii) PCT, the relevant background art 
disclosed in the document D1 is not mentioned in the description, nor is this 
document identified therein. 
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COMPENSATION OF DOPPLER SHIFT 
IN A MOBILE COMMUNICATION SYSTEM. 

The present invention relates to compensation of Doppler shift 
in a mobile communication station. 

In a mobile communication system, signals which are transmitted 
from mobile stations moving relative to a base station are 
subject to the well-known Doppler effect, which causes a 
frequency shift in the frequency received at the base station 
relative to that which was transmitted by the mobile station. 
This frequency shift is referred to herein as the Doppler shift . 
The Doppler shift is dependent upon the speed and direction of 
movement of the mobile station relative to the base station. 
Thus, the Doppler effect can provide an increase or a decrease 
in the frequency, depending on the direction of movement of the 
mobile station relative to the base station. The magnitude of 
the Doppler shift is dependent on the speed with which the mobile 
station is moving relative to the base station 1 . 

Existing mobile communication installations provide a form of 
Doppler compensation, in that the frequency detection circuitry 
within the base station which selects a particular signal on a 
particular channel can take into account a certain amount of 
Doppler shift in the signal. 

AU 664626 relates to a method and circuit arrangement for 
compensating for the Doppler shift in a radio signal propagating 
between a base station and a mobile station when the mobile 
station approaches and moves past the base station. As a mobile 
station approaches the base station with a decreasing propagating 
time delay, the propagating time delay is integrated at intervals 
to determine how the propagating time delay varies with time. 
This information is used to effect a change in the radio signal 
frequency at a particular time to compensate for the sudden 
Doppler shift as the mobile station moves past the base station. 
In a TDMA system, the interrogation occurs at intervals equal to 
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an integral number of time frames. In the GSM standard a time 
frame comprises eight consecutive time slots and a single 
transmission burst passes between a particular mobile station and 
base station in any one time frame. 

Thus, in this method of Doppler compensation decisions are made 
in response to past and incoming signals to improve the reception 
at a future time, i.e. it is a reactive system. It would be 
desirable to implement a system which can actively compensate for 
Doppler shifts in an incoming signal in real time. 

According to one aspect of the invention there is provided a 
method for compensating for Doppler shift in a signal transmitted 
between a mobile station and a base station in a mobile 
communication system, the method comprising: 

determining a channel impulse response for the channel on 
which the signal is received; . 

using the channel impulse response to estimate data bits of 
a selected portion of the received signal; 

generating a reference vector using the channel impulse 
response and the estimated data bits; 

determining a Doppler characteristic using the selected 
portion of the received signal and the reference vector; and 

using the Doppler characteristic to provide a Doppler shift 
compensation for the received signal. 

According to another aspect of the invention there is provided 
a system for compensating for Doppler shift in a signal 
transmitted between a mobile station and .a base station in a 
mobile communication system, the system comprising: 

a channel impulse response determination circuit for 
determining a channel impulse response for the channel on which 
the signal is received; 

an estimation circuit connected to receive the received 
signal on the channel impulse response and to estimate data bits 
of a selected portion of the received signal using the channel 
impulse response; 
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a reference generator for generating a reference vector 
using the channel impulse response and the estimated data bits; 

circuitry for determining a Doppler characteristic using the 
selected portion of the received signal and the reference vector; 
and 

a Doppler shift compensation circuit operable to use the 
Doppler characteristic to provide a Doppler shift compensation 
for the received signal. 

The invention is particularly applicable in TDMA mobile 
communication systems where the signal comprises a transmission 
burst. The selected portion is located in the transmission burst 
close to a zero phase offset point so that the effect of the 
Doppler characteristic is sufficiently small that it does not 
corrupt the transmitted bits.' 

According to the GSM standard, a TDMA transmission burst 
comprises a training sequence which is normally used to determine 
the channel impulse response for the channel on which the signal 
is received. This is done by a convolution of the received 
training sequence with a stored version of the training sequence. 

In existing systems, the channel impulse response is used to 
remove from the received signal the effects of the transmission 
channel on the signal, in particular multi-path and attenuation 
effects . 

According to the GSM standard, the "cleaned up" and filtered 
signal is demodulated to remove the IQ modulation by means of 
which the data in the signal was transmitted. Then, the 
demodulated signal can be decoded to generate hard bits. This 
can be done by a Viterbi technique. 

The Doppler shift is estimated from samples of the received 
signal. The estimate for Doppler shift is thus dependent on 
channel quality, and typically gets worse when the channel 
quality is poor and better as the channel quality improves. When 
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there is no Doppler shift and the channel conditions are near the 
sensitivity level of a receiver, the application of a Doppler 
compensation algorithm degrades the performance of the receiver. 

On the other hand, if there is a Doppler shift in poor channel 
conditions near the sensitivity level of the receiver, the 
receiver is not able to meet the reference sensitivity limits if 
a Doppler compensation algorithm is implemented. 

According to one embodiment of the present invention the method 
for compensating for Doppler shift in a signal transmitted 
between a mobile station and a base station in a mobile 
communication system, includes: 

detecting the quality of the received signal; and 
implementing a Doppler shift compensation in dependence on 
the detected signal quality. 

Thus in this embodiment Doppler shift compensation is used only 
or mainly in good enough channel conditions. This provides an 
increase in performance for the receiver. The technique is 
provided for ensuring that Doppler compensation does not degrade 
the sensitivity of the receiver. 

The step of detecting signal quality can include estimating the 
noise energy component of the signal. ' This can either be used 
itself to generate a Doppler' correction modification factor for 
controlling the Doppler shift compensation in dependence on the 
detected signal quality, or to generate a signal to noise ratio 
for the received signal, which would then be used to generate the 
Doppler correction modification factor. 

The Doppler compensation can be implemented as any appropriate 
user defined function of signal quality. For example, it could 
be a linear function or a step function. 

In the described embodiment, the adaptive method for Doppler 
correction applied to a GSM system first estimates the quality 
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of the channel and uses the resulting modification factor to 
scale a calculated phase difference between a reference signal 
and the actual received signal. 

For a better understanding of the present invention and to show 
how the same may be carried into effect, reference will now be 
made by way of example to the accompanying drawings in which: 

Figure 1 is a diagram of a signal burst in a mobile 
communication system; 

Figure 2 is a block diagram of circuitry for implementing 
modified Doppler shift compensation; and 

Figure 2a is a block diagram of a revised Doppler correction 
modification factor generator circuit. 

Figure 1 illustrates a normal burst in a mobile communication 
system according to the GSM standard. This figure represents a 
burst received at a base station. For a TDMA system according 
to the GSM standard, s mobile stations transmit bursts as modulated 
signals on frequency channels allocated by a base station 
controller. One frequency channel may support up to eight 
bursts, each burst associated with a respective call, where each 
call is allocated a time slot in which to send the burst. 
Further details of a TDMA system according to the GSM standard 
are not described herein because they are known to a person 
skilled in the art. 

The normal burst contains two packets of 58 bits (DATA) 
surrounding a training sequence (TRS) of 26 bits. Three tail 
bits (TS) are added at each' end of the normal burst. The 
training sequence (TRS) is a. predetermined sequence of bits which 
is sent by the mobile station (MS) and is known at the base 
station controller (BSC) . It is utilised at the base station 
controller to estimate the impulse response of the channel over 
which the burst is sent. The actual information which is 
transmitted is located in the data bits (DATA) of the burst. 
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As explained earlier, the environment through which a signal 
passes from a mobile station to a base station can vary 
considerably, depending, amongst other things, on the distance 
between the mobile station and the base station, and interference 
caused by buildings and other structures in the area. As a 
result, the signal strength and signal quality of the signal 
received at the base station varies widely. Moreover, for moving 
mobile stations, the signal received by the base station is 
subject to a Doppler shift which should be corrected. 

The circuit described herein provide a Doppler shift correction. 
The particular circuit described in relation to Figure 2 provides 
a correction only in situations where the channel conditions are 
good enough to give adequate signal quality received at the base 
station. Thus, a modification factor S c is generated dependent 
on signal quality, which is used to control the Doppler shift 
correction so that Doppler shift correction is applied only in 
good enough channel conditions. It should be appreciated that 
although the use of a modification factor S c in advantageous it 
is not essential to implementing a Doppler correction per se. 

Figure 2 illustrates a circuit 1 suitable for implementing a 
Doppler compensation in a GSM system. It should be understood 
that the various blocks in Figure 2, although illustrated as 
separate interconnected entities, do not necessarily represent 
separate physical entities, but are intended to represent 
diagrammatically the various steps which are carried out. The 
blocks could be implemented as circuits or a suitably programmed 
microprocessor may effect each of the functions which is 
individually assigned to the blocks. 

An antenna 20 receives signals 11 from the mobile stations. The 
antenna 20 is connected, via an interconnect 21, to RF circuitry 
22. This circuitry 22 operates on the received burst to 
downshift the frequency to the baseband frequency and to sample 
the burst to provide from the analogue signal digital sampled 
values. The output of RF circuitry 22 is a sampled burst r (in 



WO 98/34357 



PCT/EP97/00488 



7 

digital form) , sampled at the expected bit rate of the 
transmitted signal. -Figure 1 illustrates such a burst. The 
output of circuitry 22 is supplied along line 24 to a channel 
impulse response (C.I.R.) block 10, to a variance calculator 16 
to enable estimation of the quality of the communication channel 
(as described later) , to filtering and equalisation circuitry 12, 
to a phase difference calculator 36 and to transforming circuitry 
40 to enable the estimation and application of a Doppler shift 
correction to the burst r. 

The top part of Figure 2 illustrates the circuitry required for 
implementing the adaptive part of the system, to generate the 
Doppler correction modification factor S c . A memory 32 holds the 
training sequence TRSref which is the predetermined sequence of 
bits which is sent by the mobile station MS as a training 
sequence and received at the base station as TRS_received . The 
reference training sequence TRSref is supplied to a reference 
generator 14 and to the channel impulse response (C.I.R.) block 
10. The reference generator 14 also receives the estimated 
channel impulse response h from' the ■ channel impulse response 
block 10. 

The C.I.R. block 10 receives the burst r, including the received 
training sequence TRS_received and calculates an estimated 
channel impulse response h by calculating the cross correlation 
between the received training sequence TRSreceived and the known 
training sequence TRSref. So, 

h = xcorr (TRS_received, TRSref) (equation 1) 

It will be appreciated that, prior to effecting the cross 
correlation, the known training sequence TRSref, which is stored 
in digital form, is i,q modulated similarly to the manner in 
which the training sequence has been modulated at the MS for 
transmission, according to the GSM standard. The cross 
correlation is done in a known manner to produce a channel 
impulse response in the form of five tap values (h (i) i-f0to4 ) 
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As is known, the estimated impulse response h is used to 
calculate the expected estimate of the data in the received 
burst, as though the data- has been subject to the same average 
noise . 

The C.I.R. block also generates timing advance information r, 
which is used to determine where in the allocated time slot the 
received burst r is located. 

For each burst, the estimated channel impulse response h for that 
burst is calculated by the CIR block 10 and is supplied to 
f iltering/equalising ' circuitry which allows the data, DATA ( r ) , 
in that burst to be recovered. As is known, the 

filtering/equalising circuit 12 receives the channel impulse 
response h and timing information r for the received burst to 
allow the signal to be demodulated, filtered and decoded, to 
recover the data in a known manner. 

The reference generator 14 produces a reference vector, reffi, 
which is calculated using the convolution of the impulse response 
and the known training sequence. Thus, the reference generator 
14 performs the following calculation: 

reffi = h*TRSref (equation 2) 

In more detail, (where reffi k represents the kth sample of the 
signal reffi) 

N-l 

reffi k = Y, hi (l-2.TRS k _i) (equation 3) 

i = 0 

in which N represents the number of tap values in the estimated 
impulse response h (N = 5 in the described embodiment), and k 
runs from N-l to 25. 

The vector reffi is supplied from the reference generator to the 
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variance calculator 16. As described above, the variance 
calculator also receives the burst r, including the received 
training sequence. The variance calculator calculates a variance 
var (a 2 ) according to the following equation: 
25 

var = (£{|r k - reffi k |*2)) (equation 4) 

k=4 



ref f i_length 

The term ref f i__length is a constant representing the length of 
the reference signal, reffi. This is calculated by multiplying 
the number of samples (22) by the bit separation; 

In equation 4, the values of r k are the sampled values of the 
received training sequence for the burst r. 

It will be appreciated that each actual received sample r k will 
have a noise level which is different to the averaged estimated 
noise level derived from the channel impulse response and 
reflected in the reference samples reffi k . Thus, the variance 
gives an indication of the level of noise energy actually 
received, and thus signal quality. 

The output a 2 of the variance calculator 16 is supplied to a 
Doppler correction modification factor circuit 18. The Doppler 
correction modification factor circuit 18 uses the calculated 
variance o 2 to generate a modification factor S c using a function 
which can be determined by a user. In the embodiment shown in 
Figure 2, the Doppler correction modification factor circuit 
generates the modification S c as a function of the variance a 2 , 
for example a linear function or a non-linear function such as 
a step function. 

In another embodiment of the invention the value of S c is 
calculated in dependence upon a signal to noise ratio (SNR) of 
the channel. Figure 2a illustrates circuitry for a Doppler 
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correction modification circuit 18' for implementing such an 
embodiment. A SNR calculator 42 receives the calculated variance 
a 2 from the variance calculator 16 and an energy value E. The 
signal energy E may be determined from the calculated signal 
reffi in accordance with the following: 
25 

E = £ |reffi k | 2 (equation 4a) 

k=4 

ref f i_length 

Alternatively the tap values can be used to derive the signal 
energy, E, of the channel according to the following: 

4 

E = E(h(i) ) 2 
i = 0 



The SNR value is calculated by the SNR calculator 42 in 
accordance with the following equation: 

SNR = E 



v 2 (equation 4b) 

The above referenced technique is described in our earlier 
Application No. PCT/FI96/00461 , the contents of which are herein 
incorporated by reference. 

The SNR value is then supplied to a revised modification factor 
generating circuit 44. This circuit 44 calculates the value of 
the modification factor S c as a function of the SNR value, for 
example a linear function or a : non-linear function such as a step 
function. 

The lower part of the circuit in Figure 2 illustrates in block 
diagram form a system for implementing a Doppler shift 
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correction. However, it will be apparent that ■ the Doppler 
correction adaptive part described above could be used with other 
implementations of Doppler correction. 

The main function of the adaptive part of the circuit described 
above is to provide a Doppler correction modification factor S c 
which is based on the channel conditions, so that Doppler 
correction is used only or mainly in good enough channel 
conditions. This is useful in any type of Doppler correction. 
A particular Doppler correction circuit is discussed in the 
following. 

An equalisation circuit 30, for example a Viterbi equaliser, 
receives the filtered, demodulated and equalised signal DATA ( r ) 
from the. filtering and equalisation circuitry 12. The 
equalisation circuit 30 operates on a part of the data sequence 
DATA of the burst (that part having been derived from ESTIM . BLOCK 
in Figure 1) to estimate and output bits which were sent from the 
mobile station MS. This output is referred to herein as 
estim_bits, and they run from k=j to k=j+n. The equalisation 
circuit 30 operates to make decisions of the bits as in known 
mobile communication systems and thus it will not be described 
further herein. 

The estimated bit decisions estim bits are supplied to a 
reference circuit 34. The reference circuit 34 generates a 
reference vector ref by using a convolution of the estimated bit 
decisions and the estimated impulse response h, according to the 
following equation : 

ref = estim_bits * h (equation 5) 

Thus, the reference vector ref comprises a set of samples ref k 
(k=j->j+n) , each having real and imaginary values. The reference 
vector ref is supplied to a phase difference calculator 36. As 
described earlier, the phase difference calculator 36 also 
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receives the received burst r. As known in the art, the received 
burst comprises samples r k each having real and imaginary values. 

The phase difference calculator uses a value t_dif f which 
represents the time between a zero phase offset point h_time and 
the -middle of the estimation block, as illustrated in Figure 1. 
The zero phase offset point h_time is a zero phase offset point 
inside the training sequence where the calculated impulse 
response is true. In practice, this is typically the middle of 
the training sequence. The value of t_diff and. h_time can be 
determined during the design phase of the system., and a constant 
thereafter. Of course, they could be reprogrammed if necessary 
during use of the system. 

In addition, the location (j) of the beginning of the estimation 
block estim__block and its length (n) is determined and programmed 
into the equalisation circuit 30. The estimation block is 
selected so that the Doppler offset has not yet corrupted the 
received bits. 

The phase change per bit duration (ph_diff) resulting from the 
Doppler effect (the Doppler characteristic) is calculated from 
the reference signal ref and the actual received signal r by the 
phase difference calculator 36 according to one of the following 
equations : 

ph_diff = 1 tan' 1 ^ imacy (r K - ref * h ) 1 

t_diff lX k real (r k -ref* k ) J (equation 6) 

ph_diff = 1 -tan" 1 jimag (r . • ref * K ) I 

t_diff k lreal(r k *ref* k ) //length(k) 

(equation 7) 

where k runs from j to j+n, and where length (k) represents the 
amount of different k values in the summation, i.e. n. 

A Doppler correction circuit 38 is then used to correct the 
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estimated Doppler shift from the received samples. The Doppler 
correction circuit 38 receives the zero phase offset point 
h_time and the Doppler correction modification factor S c . 
Furthermore, it receives the calculated phase difference from the 
phase difference calculator 36. Knowing that the point h_time 
has zero phase offset, the actual Doppler phase shift <f> can be 
calculated for each bit as follows: 

0k = S c * ph_diff • (k-h_time) (equation 8) 

where k is a bit index of the received sample r. When the index 
k<h_time the phase shift has an opposite sign to when k>h_time. 
It will be apparent that if the Doppler correction circuitry is 
operating in independence of the adaptive circuitry which 
produces Sc, Sc will default to 1 in equation 8. 

Transforming circuitry 40 is then used to implement the Doppler 
shift correction on the received burst r to produce a corrected 
signal. The transforming circuitry receives the estimated 
Doppler shift vector <j> (comprising the <f> k values) and sampled 
values of the received burst r. It performs a CORDIC operation 
to correct for the Doppler shift of each sample, according to the 
following operation. 

( N_real_sample(k) U( cos<p k sintf> k ) (real(r k )\ 

I N_imag_sample (k) / \ sin0 k cos0 k / \imag(r k )/ 

(equation 9) 

The Doppler shift corrected vector DCV which is output from the 
transforming circuitry 40 is supplied to the filtering/equalising 
circuit 12, so that the Doppler corrected signal is used to 
recover the data from the signal. 

As part of the Doppler correction technique described above, bit 
decisions for the estimation block have already been made. These 
constitute part of the data. It is thus not necessary to 
estimate the same bits again, although this could be done. 
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Instead, the equalisation circuit 30 can be stopped at the end 
of the estimation block and the current state preserved. Next, 
the Doppler correction is performed for the remaining bits and 
the Viterbi equalisation is then executed to the end of the time 
slot on the Doppler corrected bits. In the second part, the 
Doppler correction may be done to the first part of the time slot 
and the Viterbi estimation can then be performed for the data 58 
in the first part ' of the time slot. This method reduces the 
required calculations in the receiver. 

It is possible to implement a limit on the value of the phase 
difference phjiiff, so that if the phase difference is below a 
certain threshold, no correction is performed. 

A typical environment where Doppler correction could be used is 
fast trains or motorways. In that situation, it is likely that 
there would be a line of sight' path from the base station to the 
mobile station, so that the Doppler correction would be a 
constant value between different time slots if the velocity of 
the mobile station is the same. In that case, an average value 
of the phase difference could be calculated from several 
different time slots, and this average value could be used as a 
correction value. 

In the case of a diversity receiver having a plurality of 
different branches, the phase differences can be calculated for 
all branches using the same estimates as the received samples, 
with each branch using its own impulse response. 
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■ CLAIMS : 

1. A method for "compensating for Doppler shift in a signal 
transmitted between a mobile station and a base station in a 
mobile communication system, the method comprising:. 

determining a channel impulse response for the channel on 
which the signal is received; 

using the channel impulse response to estimate data bits of 
a selected portion of the received signal; 

generating a reference vector using the channel impulse 
response and the' estimated data bits; 

determining a Doppler characteristic using the selected 
portion of the received signal and the reference vector; and 

using the Doppler characteristic to provide a Doppler shift 
compensation for the received signal. 

2. A method according to claim 1, wherein the signal comprises 
a transmission burst in a TDMA mobile communication system. 

3. A method according to claim 1, wherein the selected portion 
is located in said transmission burst close to a zero phase 
offset point, whereby it has been substantially unaffected by the 
Doppler characteristic. 

4. A method according to claim 2 or 3 , wherein the channel 
impulse response is determined from a training sequence in said 
transmission burst. 

5. A method according to claim 1, 2, 3 or 4, wherein the data 
bits are. estimated by: 

using the channel impulse response to remove from the 
received signal the effects of the transmission channel for the 
signal ; 

demodulating the resulting signal; and 

decoding the selected portion of the demodulated signal to 
estimate data bits. 
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6. A method according to claim 5, wherein the decoding step is 
done by a Viterbi method.- 

7. A method according to any preceding claim, wherein the step 
of using the Doppler characteristic comprises determining for 
each of a plurality of indexed samples of the received signal the 
Doppler phase shift for that sample based on the Doppler 
characteristic and the location of the induced sample within the 
received signal. 

8. A method according to claim 7, where the samples are indexed 
at the bit rate of the received signal. 

9. A system for compensating for Doppler shift in a signal 
transmitted between a mobile station and a base station in a 
mobile communication system, the system comprising: 

a channel impulse response determination circuit for 
determining a channel impulse response for the channel on which 
the signal , is received; 

an estimation circuit connected to receive the received 
signal on the channel impulse response and to estimate data bits 
of a selected portion of the received signal using the channel 
impulse response; 

a reference generator for generating a reference vector 
using the channel impulse response and the estimated data bits; 

circuitry for determining a Doppler characteristic using the 
selected portion of the received signal and the reference vector; 
and 

a Doppler shift compensation circuit operable to use the 
Doppler characteristic to provide a Doppler shift compensation 
for the received signal. 

10. A system according to claim 9, wherein the estimation 
circuit comprises filtering an equalisation circuit operable to 
remove from the received signal the effects of the transmission 
channel for the signal by using the channel impulse response; 

demodulation circuitry for demodulating the resulting 
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signal; and 

decoding circuitry for decoding the selected portion of the 
demodulated signal to estimate the data bits. 

11. A system according to claim 10, wherein the decoding 
circuitry is a Viterbi decoder. 
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AMENDED CLAIMS 

[received by the International Bureau on 3 March 1998 (3.03.98); 
original claims 1,2,3,9 and 10 amended; remaining 
claims unchanged (3 pages)] 

1. A method for compensating for Doppler shift in a signal 
transmitted between a mobile station and a base station in a 
mobile communication system, said signal comprising a sequence 
of transmission bursts, the method comprising: 

determining for each transmission burst a channel impulse 
response for the channel on which the signal is received; 

using the channel impulse response to estimate data bits of 
a selected portion of the transmission burst; said selected 
portion being located in said transmission burst close to a zero 
phase offset point; 

generating a reference vector using the channel impulse 
response and the estimated data bits; 

determining a Doppler characteristic in the form of a phase 
shift per bit duration using the selected portion of the received 
signal and the reference vector; and 

using the Doppler characteristic determined for the selected 
portion to provide a Doppler shift compensation- . for the 
transmission burst. 

2. A method according to claim l, when used in a TDMA mobile 
communication system. 

3. A method according to claim 1, wherein the Doppler shift 
compensation <f> k for each bit is calculated as follows, wherein k 
is the bit index, h_time is the point of zero phase offset and 
ph_diff is the phase shift per bit duration, 

ph__diff • (k-h_time) . 

4. A method according to claim 2 or 3 , wherein the channel 
impulse response is determined from a training sequence in said 
transmission burst. 

5. A method according to claim 1, 2/ 3 or 4, wherein the data 
bits are estimated by: 
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using the channel impulse response to remove from the 
received signal the effects of the transmission channel for the 
signal; 

demodulating the resulting signal; and. 

decoding the selected portion of the demodulated signal to 
estimate data bits. 

6. A method according to claim 5, wherein the decoding step is 
done by a Viterbi method. 

7. A method according to any preceding claim, wherein the step 
of using the Doppler characteristic comprises determining for 
each of a plurality of .indexed samples of the received signal the 
Doppler phase shift for that sample based on the Doppler 
characteristic and the location of the induced sample within the 
received signal. 

8. A method according to claim 7 , where the samples are indexed 
at the bit rate of the received signal. 

9. A system for compensating for Doppler shift in a signal 
transmitted between a mobile station -and a base station in a 
mobile communication system, said signal comprising a sequence 
of transmission bursts, the system comprising: 

a channel impulse response determination circuit for 
determining for each transmission burst a channel impulse 
response for the channel on which the signal is received; 

. an estimation circuit connected to receive the received 
signal on the channel impulse response and to estimate data bits 
of a selected portion of the transmitted burst using the channel 
impulse response, said selected portion being located in said 
transmission burst close to a 2ero phase offset point; 

a reference generator for generating a reference vector 
using the channel impulse response and the estimated data bits; 

circuitry for determining a Doppler characteristic in the 
form of a phase change per bit duration using the selected 
portion of the received signal and the reference vector; and 
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a Doppler shift compensation circuit operable to use the 
Doppler characteristic to provide a Doppler 1 shift compensation 
for the transmission burst. 

10. A system according to claim 9, wherein the estimation 
circuit comprises: 

filtering and equalisation circuitry operable to remove from 
the received signal the effects of the transmission channel for 
the signal by uSing the channel impulse response; 

demodulation circuitry for demodulating- the resulting 
signal; and 

decoding circuitry for decoding the selected portion of the 
demodulated signal to estimate the data bits. 

■11. A system according to claim 10, wherein the decoding 
circuitry is a Viterbi decoder. 
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Die Erfindung betrifft ein Verfahren nach dem Oberbegriff des Patentanspruches 1 . 
Die Erfindung findet Verwendung bei digitalen Kommunikationsnetzen, insbesondere Mobilfunksystemen. 
Das Verfahren ist etwa fur TDMA (Time Division Multiplexing Access) Systeme, z.B. das GSM (Group 
Speciale Mobile) System geeignet. 
5 Im Mobilfunk ist die Kanalimpulsantwort durch den Dopplereffekt zeitlich veranderlich. Es stellt sich das 
Problem, diese Anderungen in einem Empfanger zu schatzen und zu berucksichtigen. 
Vorbekannte Losungen sind beispielsweise: 

a) die Kanalimpulsantwort wird im gesamten Zeitschlitz als konstant angenommen. 

b) Die Kanalimpulsantwort wird vollstandig im Zeitschlitz nachgefuhrt. Dies bedeutet, die zeitliche 
io Veranderung der Mehrwegeausbreitung und damit die Form der Kanalimpulsantwort fur alle Zeitpunkte 

im Zeitschlitz zu schatzen. Die Nachfuhrung kann mit Hilfe eines Kalmanfilters geschehen. In jedem 
Falle ist entweder eine Entscheidungsruckkopplung erforderlich oder es sind weitere Stutzstellen (d.h. 
dem Empfanger bekannte feste Sequenzen) im Zeitschlitz erforderlich. Auch ist eine Interpolation der 
Kanalimpulsantwort im ganzen Zeitschlitz moglich, falls hintereinanderliegende Zeitschlitze auf der 
75 gleichen Frequenz entsprechend dicht liegen und die Geschwindigkeit nicht zu gro/3 ist (Literatur: Hoher, 
P.: Koharenter Empfang trelliscodierter PSK-Signale auf frequenzselektiven Mobilfunkkanalen - Entzerrer, 
Decodierung und Kanalparameterschatzung, Fortschritt-Bericht VDI, Dissertation, Universitat Kaiserslau- 
tern, Mai 1990). 

Die bekannten Verfahren haben den Nachteil, da/J bei hohen Geschwindigkeiten sich eine entsprechend 
so hohe Bitfehlerrate ergibt und dafl entsprechende Verfahren sehr aufwendig sind. 

Im Falle von Entscheidungsruckfuhrung reagieren diese Verfahren sehr empfindlich auf Bitfehler. Bei 
kleiner Geschwindigkeit und hohen Storpegeln ergibt sich eine grofle Degradation der Empfindlichkeit im 
Vergleich zu einem einfachen Empfanger, bei dem die Kanalimpulsantwort im Zeitschlitz als konstant 
angenommen wird. 

25 Werden Stutzstellen im Zeitschlitz verwendet, so bedeutet dies naturlich eine Verkleinerung der 
Ubertragungskapazitat. Die Benutzung von Nachbarzeitschlitzen verbietet die Benutzung eines Frequenz- 
sprungverfahrens, welches gro/te Vorteile bei Fading aufweiflt und die TDMA Struktur mu/3 an die maximale 
Geschwindigkeit angepasst werden. 

Der Erfindung liegt deshalb die Aufgabe zugrunde, ein Verfahren anzugeben, bei dem eine Dopplerkor- 

30 rektur und Tragerfrequenzkorrektur bei Ubertragung in Zeitmultiplex-Systemen mit linearen, digitalen 
Modulationen mit konstanter einhullender unter Minimierung des Aufwandes bei gleichzeitig hoher Emp- 
fangsqualitat erreicht wird. Diese Aufgabe wird gelost durch die im kennzeichnenden Teil des Patentan- 
spruch 1 angegebenen Merkmale. Vorteilhafte Ausgestaltungen und/oder Weiterbildungen sind den Unter- 
anspruchen zu entnehmen. 

35 Es wird davon ausgegangen, da/3 die Zeitschlitze des TDMA-Systems in der Mitte eine dem Empfanger 
bekannte Trainingssequenz enthalten, mit deren Hilfe die Kanalimpulsantwort gemessen werden kann. Die 
Erfindung besteht zunachst aus der Vorabdemodulation des mittleren Teiles eines Zeitschlitzes. Danach 
kann eine mittlere Phasendrehung berechnet werden. Diese beschreibt in guter Naherung auch bei 
Mehrwegeausbreitung und bei beliebigem Dopplerspektrum (z.B. symmetrisches Dopplerspektrum) die 

40 Anderung der Kanalimpulsantwort innerhalb eines Zeitschlitzes. Die endgultige Demodulation erfolgt dann in 
mehreren Bereichen innerhalb eines Zeitschlitzes. In jedem Bereich wird hierzu die Kanalimpulsantwort als 
konstant angenommen und entsprechend der mittleren Phasendrehung berechnet. Nach Mittelung uber 
mehrere Zeitschlitze kann die Phasendrehung auch fur die Berechriung der Stellgro/te fur die analogen 
Tragerfrequenzerzeugung verwendet werden. 

45 Die Erfindung hat den Vorteil, dafl lediglich kleine Bitfehlerraten bei hohen Geschwindigkeiten auch bei 
Mehrwegeausbreitung mit einer maximalen Mehrwegeverzogerurig, welche gro/Jer als die Symboldauer ist, 
auftreten. Gleichzeitig ist nur ein geringer Mehraufwand erforderlich im Vergleich zu einem einfachen 
Verfahren, welches von einer konstanten Kanalimpulsantwort uber den gesamten Zeitschlitz ausgeht Durch 
die Dopplerkorrektur ergibt sich keine Verschlechterung der Empfindlichkeit bei kleinen Geschwindigkeiten. 

so Weitere Stutzstellen auflerhalb der Trainingssequenz werden nicht benotigt. Bei Intersymbolinterferenz kann 
wie bei einem einfachen Verfahren z.B. ein Viterbi Entzerrer verwendet werden, der von einer konstanten 
Autokorrelationsfunktion der Kanalimpulsantwort uber den gesamten Zeitschlitz ausgeht. Die Erfindung 
benotigt kein Gedachtnis uber mehr als einen Zeitschlitz und damit ist ein Frequenzsprungverfahren 
moglich. 

55 Die Erfindung wird im folgenden anhand von AusfGhrungsbeispielen beschrieben unter Bezugnahme auf 
schematische Zeichnungen. 

Es wird ein digitales Modulationsverfahren (z.B. Phase Shift Keying PSK oder Minimum Shift Keying 
MSK) und ein linearer Ubertragungskanal betrachtet. Es ist allgemein bekannt, da£ das Empfangssignal y(t) 
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dann durch die Faltung einer Kanalimpulsantwort p(t,r) mit den gesendeten komplexen oder reellen 
Symbolen a« beschrieben werden kann: 

y{t) = / pit, t) Y\ &k6{t - r - kT)dr 



Die Kanalimpulsantwort gibt sowohl die Eigenschaften der Modulation wie des .Funkkanales wieder 

70 (Mehrwegeausbreitung). Sie beschreibt die Antwort der Ubertragungsstrecke zum Zeitpunkt t auf einen 
DiraoSto/3 als Eingangssignal in Abhangigkeit der Verzogerung r. Die Kanalimpulsantwort wird zunachst als 
zeitveranderlich angesehen. Durch Do'pplerverschiebungen andert sich z.B. im Mobilfunk die Form der 
Kanalimpulsantwort nach einer gewissen Wegstrecke vollstandig. 1st die maximale Mehrwegeverzogerung 
klein gegenuber der Symboldauer, so kann die zeitliche Veranderung der Kanalimpulsantwort durch eine 

75 Amplituden- und Phasenanderung beschrieben werden. Die Form der Kanalimpulsantwort hangt dann nicht 
mehr vom Funkkanal ab. Es wird nun ein TDMA System betrachtet, bei dem fur die Mehrwegeverzogerung 
mehrere Symboflangen zugelassen werden, Jeder Zeitschiitz besitzt in der Mitte eine teste dem Empfanger 
bekannte Symbolsequenz (Trainingssequenz). Mit ihr wird die Kanalimpulsantwort in der Mitte des Zeit- 
schlitzes gemessen werden. Fur eine richtige Demodulation muJ3 jedoch auch die zeitliche Anderung der 

20 Kanalimpulsantwort innerhalb eines Zeitschlitzes bekannt sein. Auch dann, wenn die Mehrwegeverzogerung 
die Dauer eines Symboles deutlich Uberschreitet, wird die zeitliche Anderung der Kanalimpulsantwort immer 
noch durch eine mittlere Amplituden- und Phasenanderung beschrieben . Voraussetzung ist, da/3 der 
Zeitschiitz nicht zu lang ist, damit die Anderung der Kanalimpulsantwort nicht zu gro/3 wird. Mit dieser 
Voraussetzung ist der Fehler der obigen Naherung gering. Es gibt also einen Wertebereich fur die Lange 

25 eines Zeitschlitzes in dem die Annahme einer zeitinvarianten Kanalimpulsantwort zu merklichen Verlusten 
der Leistungsfahigkeit bei der Demodulation fuhrt. Die Annahme einer mittleren konstanten Amplituden- und 
Phasenanderung der Kanalimpulsantwort innerhalb eines Zeitschlitzes fuhrt dann zu einer deutlichen 
Verbesserung. 

Bei Phasenmodulationsarten mit konstanter Einhullender steckt die Information nur in der Phase und auch 
30 bei einer Anderung der Amplitude der Kanalimpulsantwort kann ohne Kenntnis dieser Amplitudenanderung 
noch richtig demoduliert werden. Deswegen reicht es dann hier aus, nur die mittlere Anderung der Phase 
zu betrachten. 

Es wird angenommen, dafl das Empfangssignal bandbegrenzt ist und daher abgetastet werden kann. 
Die Abtastrate soli n mal die Symbolrate betragen. Das abgetastete Empfangssignal wird nun durch die 
35 diskrete Faltung der Symbole mit der abgetasteten Kanalimpulsantwort beschrieben werden. 

Im+k = 1^-02^+* * = °. * ' ' , n - 1 

40 

Y ni + k sind die empfangenen komplexen Abtastwerte, P t , n i + k ist die Kanalimpulsantwort und die ai sind die 
gesendeten komplexen (z.B. 4,8 PSK) oder reellen Symbole (z.B. 2PSK oder MSK). 

Zunachst wird die Kanalimpulsantwort in der Mitte des Zeitschlitzes mit Hilfe einer festen dem 
45 Empfanger bekannten Sequenz (Trainingssequenz) geschatzt. Dann kann in einem mittleren Teil des 
Zeitschlitzes (siehe FIG. 1) mit Hilfe eines digitalen Matched Filters (an die gemessene Kanalimpulsantwort 
angepaGt) eine Symbolschatzung vorgenommen werden (Vorabdemodulation). 

k 



Die empfangenen Symbole lassen sich durch eine Harddecisio'n des Ausgangs des Matched Filters mfoi 
55 bestimmen. Es werden nun rechts und links der Trainingssequenz zwei Bereiche definiert, deren Symbole 
nun bekannt sind. Diese Bereiche werden im folgenden Ersatzpraambeln genannt. Es wird nun jeweils die 
mittlere Phasendrehung der Kanalimpulsantwort zwischen dem Bereich einer Ersatzpraambel und dem 
Bereich der Trainingssequenz geschatzt. Durch Mittelung beider Winkel wird ein endgiiltiger Wert fur die 
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Phasendrehung bestimmt. 

Im folgenden wird nun die Schatzung der Phasendrehung mit Hilfe einer Ersatzpraambel und der 
bereits bekannten Kanalimpulsantwort beschrieben. Es wird vorausgesetzt, da/3 die Ersatzpraambel so kurz 
sein soil, da/3 man die Phase <f> innerhalb der Ersatzpraambel als konstant ansetzen kann. Das Empfangssi- 
5 gnal ist dann durch folgende Gleichung gegeben. 
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Nach dem Maximum Likelihood Prinzip kann jetzt die Phase geschatzt werden. Der folgende Ausdruck mu/3 
dazu minimiert werden (der Rauschanteil aufeinanderfolgender Abtastwerte wird als unkorreliert angesehen, 
wei/Jes Rauschen): 



n-1 



i k=zQ I 



Die Minimierung erfolgt bezuglich Alle anderen Grotfen sind bereits bekannt und damit test Die 
Minimierung des obigen Ausdrucks laflt sich dann durch die Maximierung des folgenden Ausdrucks 
erreichen: 



*e{e-^££p; (+fc £y ni+fc a:_J 

i k=o > i 



35 Der Realteil wird maximal, falls fUr $ der Phasenwinkel der Summe gewahlt wird. Dies bedeutet, da/3 zur 
Berechnung der Phase J die konjugiert komplexen vorabdemodulierten Symbole der Ersatzpraambel mit 
dem Empfangsabtastwerten korreliert werden mussen. 

Dies entspricht einer erneuten Schatzung der Kanalimpulsantwort. Die bereits bekannte Kanalimpulsant- 
wort wird nun mit der neu geschatzten konjugiert komplexen Kanalimpulsantwort Abtastwert fur Abtastwert 

40 multipliziert und aufakkumuliert. 

Da entsprechend den obigen Voraussetzungen die Phase Goer einen gro/Jeren Bereich als konstant 
angesehen werden kann, braucht auch bei der endgiiltigen Bestimmung der empfangenen Symbole der 
Phasenwinkel nicht Symbol fur Symbol nachgefuhrt werden. Es reicht aus, den Zeitschlitz in mehrere 
Bereiche einzuteilen, in denen die Phase jeweils wieder als konstant angesehen werden kann (siehe FIG. 1). 

45 Die Symbole werden wieder mit Hilfe eines Matched Filters geschatzt. Die Kanalimpulsantwort wird fur 
jeden Bereich entsprechend der geschatzten Drehung und dem Abstand von der Trainingssequenz 
bestimmt: 



A ist der Abstand einer Ersatzpraambel zur Trainingssequenz und d ist der mittlere Abstand eines 
55 Bereiches fur das endgiiltige Matched Filter zur Trainingssequenz. Zur Korrektur der Intersymbolinterferenz 
. konnen bekannte Standardverfahren angewendet werden (z.B. Viterbi Entzerrer). Es mu/3 hierbei nicht mehr 
zwischen verschiedenen Bereichen im Zeitschlitz unterschieden werden, da die Autokorrelationsfunktion der 
Kanalimpulsantwort nicht von einer Phasendrehung abhangig ist. 
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FIG. 2 zeigt ein Btockschaltbild des vollstandigen Empfangers. 

In jedem Empfanger besteht das Problem die Tragerfrequenz analog nachzufuhren. Dies kann mit dem 
obigen Verfahren kombiniert werden. Die geschatzte Phasendrehung gibt nicht nur die Doppierverschiebung 
wieder, sondern nach entsprechender Mittelung auch den Versatz der Tragerfrequenz (siehe FIG. 3). 

Im folgenden werden einige Simulationsergebnisse fur einen GSM Empfanger nach obigem Verfahren 
gegeben. Das endgQItige Matched Filter arbeitet hierbei in drei verschiedenen Bereichen und als Modula- 
tion wird GMSK verwendet (binare Modulation). Zum Vergleich werden Ergebnisse fur einen Empfanger 
dargestellt, der eine konstante Kanalimpulsantwort im gesamten Zeitschlitz annimmt. Als Ausbreitungsmo- 
dell wird das sogenannte "Rural Area" Profil nach den COST 207 Spezifikationen verwendet. Dieses Profil 
ist im GSM System eines der Referenzprofile zur Spezifikation der erlaubten Bitfehlerraten. Angegeben wird 
im folgenden die Restfehlerwahrscheinlichkeit des Demodulators bei Empfang eines Signales mit dem 
obigen Mehrwegeprofil ohne weitere Storungen durch Rauschen in Abhangigkeit der Geschwindigkeit V. 



V in km/h 


50 




250 




500 


Neues Verfahren 


3.2 • 


10"* 


5.6 • 


10" 3 


2.1 * 10" 2 


Verfahren mit konst. angenommener Kanalimpulsantwort 


3.8 • 


10 -4 


9.0 • 


10" 3 


3.3 • 10" 2 
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Die Verbesserung durch das neue Verfahren insbesondere bei hohen Geschwindigkeiten ist deutlich zu 
20 sehen. 

Die folgende Tabeile gibt die Bitfehlerrate des Demodulators bei Empfang des obigen Profiles bei 50 
km/h und mit Rauschen (E b /N 0 = 6 dB) an: 



25 



Neues Verfahren 

Verfahren mit konst. angenommener Kanalimpulsantwort 



6.8 • 10" 2 
6.8 • 10" 2 
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Daraus geht deutlich hervor, da/3 sich durch die Dopplerkorrektur bei kleinen Geschwindigkeiten keine 
Verschlechterung der Empfindlichkeit ergibt. 

Patentahspruche 



1. Zeitmultiplex- Verfahren zur Bestimmung der mittleren Phasenanderung des Empfangssignals, bei dem 
fur die Mehrwegeyerzogerung mehrere Symbollangen verwendet werden und jeder Zeitschlitz in der 

35 Mitte eine feste dem Empfanger bekannte Symbolsequenz, sog. Trainingssequenz, besitzt dadurch 
gekennzeichnet, 

- daj3 die Kanalimpulsantwort innerhalb eines Zeitschlitzes mit einer festen dem Empfanger 
bekannten Trainingssequenz geschatzt wird, 

- da/3 im mittleren Teil des Zeitschlitzes eine an die geschatzte Kanalimpulsantwort angepaflte 
40 Symbolschatzung durch Demodulation vorgenommen wird, 

- da/J mit den empfangenen Symbolen rechts und links der Trainingssequenz zwei Bereiche 
definiert werden, sog. Ersatzpraambeln, 

- da/3 zwischen dem Bereich der Ersatzpraambeln und dem Bereich der Trainingssequenz die 
mittleren Phasendrehungen der Kanalimpulsantwort geschatzt werden und durch Mittelung Ober 

45 die Winkel der Schatzwert fur die genaue Phasendrehung verbessert wird. 

2. Verfahren nach Anspruch 1 , dadurch gekennzeichnet, 

- da/3 die Phasendrehung der Kanalimpulsantwort mit einem Empfangssignal mit konstanter Phase 
innerhalb der Ersatzpraambel dadurch bestimmt wird, da/3 die konjugiert komplexen, vorab 

so demodulierten Symbole der Ersatzpraambel mit den Empfangssignalabtastwerten korreliert wer- 

den und die Kanalimpulsantwort erneut geschatzt wird, und 

- da/3 die bekannte Kanalimpulsantwort mit der neu geschatzten, konjugiert komplexen Kanalim- 
pulsantwort Abtastwert fOr Abtastwert multipliziert und aufakkumuliert sind. 



55 3. Verfahren nach Anspruch 2 t dadurch gekennzeichnet, da/3 der Zeitschlitz in mehrere kurze Bereiche 
annahernd konstanter Phase eingeteilt wird, 

- da/J die Kanalimpulsantwort fClr jeden Bereich entsprechend der geschatzten Phasendrehung und 
dem Abstand von der Trainingssequenz durch Demodulation bestimmt wird und iiber die mittlere 
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Phasenanderung die Anderung der Kanalimpulsantwort ermittelt wird. 

Verfahren nach einem der vorherigen Anspruche, dadurch gekennzeichnet, daJ3 mit der geschatzten 
Phasendrehung durch entsprechende Mittelung der Versatz der Tragerfrequenz bestimmt wird und 
Liber einen Regelkreis die Tragerfrequenz reguliert wird. 
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© Zeitmultiplex-Verfahren zur Bestimmung der mittleren Phasenanderung eines Empfangssignals. 



© Die Erfindung betrifft ein Zeitmultiplex-Verfahren, 
mit dem die mittlere Phasenanderung eines Emp- 
fangssignals bestimmt wird. Im Mobilfunk ist die 
Kanalimpulsantwort durch den Dopplereffekt zeitlich 
veranderlich. Es stent sich das Problem, diese Ande- 
rungen in einem Empfanger zu schatzen und zu 



berucksichtigen. Es wird ein einfaches Verfahren an- 
gegeben, welches die Anderungen der Kanalimpuls- 
antwort durch eine mittlere Phasenanderung be- 
schreibt. Das neue Verfahren ist fur TDMA Systeme 
wie zum Beispiel das GSM System geeignet. 
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(54) Compensation of doppler errors, particular in cellular networks 



(57) In accordance with one embodiment of the in- 
vention, a method of compensating for Doppler error in 
a wireless communications system employing Viterbi 
decoding comprises the steps of: for each signal sample 
in a first predetermined-sized grouping of received sig- 
nal samples, performing a parallel Viterbi update and 
short symbol decode; and for a second predetermined- 
sized grouping, forming by pipeline processing an esti- 
mate of the Doppler error in accordance with the parallel 
short traceback decoding performed for the first group- 
ing, and adjusting each signal sample in the second 



grouping in accordance with the estimated Doppler er- 
ror. 

In accordance with another embodiment of the in- 
vention, a Viterbi traceback reconstructed signal sample 
index comprises: a state counter, a traceback shift reg- 
ister (TBSR); a signal reconstruction table; and a com- 
parator coupled in a configuration so as to provide the 
sign bit to the traceback shift register from a comparison 
of binary digital signals. The state counter is coupled so 
as to provide digital signals to the traceback shift register 
and the traceback shift register is coupled so as to pro- 
vide digital signals to the signal reconstruction table. 
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Description 

Technical Field 

The present invention relates to communications 
and, more particularly, to wireless communications sys- 
tems. 

Background of the Invention ' 

Several standards are being employed with respect 
to signaling standards for digital cellular telephony 
worldwide. One such standard is Europe's global sys- 
tem for mobile communications (GSM), described, for 
example, in the ETSI/GSM Specifications. Likewise, an- 
other standard, employed in the United States, is IS54. 
These signaling standards are typically intricate. Typi- 
cally, a signal burst is transmitted via a wireless medium 
comprising a plurality of differentially encoded symbols, 
such as bits or binary digital signals. The encoded sym- 
bols are typically transmitted using phase shift keying in 
the baseband, such as, for example, minimum phase 
shift keying (MSK). This baseband signal form of the sig- 
nal burst is typically then transmitted via a wireless me- 
dium using a radio frequency (RF) carrier. 

In wireless communications, such as digital cellular 
telephony, transmission and reception typically occurs 
between stations, at least one of which is in motion. Fur- 
thermore, the velocity of the mobile station may change 
over time. As is well-known, such relative motion be- 
tween transmitting and receiving stations, such as be- 
tween a mobile station and a base station, may result in 
a Doppler shift of the frequency of the signal being trans- 
mitted. This Doppler shift may therefore result in a phase 
or frequency error in the received signal. Because of 
such a Doppler shift, the integrity of the signal being 
transmitted may be corrupted at the receiving end of the 
communications system. 

Although various processes for compensating for 
such errors due to a Doppler shift are known, one prob- 
lem associated with conventional approaches is the im- 
pact upon available signal processing resources at the 
receiving station. Typically, a receiving station, such as 
a mobile station, has limited signal processing capabil- 
ity. Therefore, exhaustive approaches to performing 
Doppler calculations and signal correction may exceed 
or at least bottleneck available resources when such re- 
sources are needed to continually process additional 
signals. A need therefore exists for a method of com- 
pensating for Doppler error while reducing associated 
bottlenecks for the available signal processing resourc- 
es of a receiving station. 

Summary of the Invention 

In accordance with one embodiment of the inven- 
tion, a method of compensating for Doppler error in a 
wireless communications system employing vlterbi de- 



coding comprises the steps of: for each signal sample 
in a first predetermined-sized grouping of received sig- 
nal samples, performing a parallel vlterbi update and 
short symbol decode; and for a second predetermined- 

5 sized grouping, forming by pipeline processing an esti- 
mate of the Doppler error in accordance with the parallel 
short traceback decoding performed for the first group- 
ing and adjusting each signal sample in the second 
* grouping in accordance with the estimated Doppler er- 

10 ror. 

In accordance with another embodiment of the in- 
vention, a Viterbi traceback reconstructed signal sample 
index comprises: a state counter, a traceback shift reg- 
ister (TBSR); a signal reconstruction table; and a corn- 
's parator coupled in a configuration. so as to provide the 
sign bit from a comparison of binary digital signals to the 
TBSR. The state counter is coupled so as to provide dig- 
ital signals to the TBSR and the TBSR is coupled so as 
to provide digital signals to the signal reconstruction ta- 
20 ble. 

Brief Description of the Drawings 

The invention will be described by way of example 
25 with reference to the accompanying drawings in which: 
FIG. 1 is a schematic diagram of one embodiment 
of a system employing a method of compensating for 
Doppler error in a wireless communications system in 
accordance with the invention. 
30 FIG. 2 is plot of complex signals in the Inphase- 
Quadrature (l-Q) plane illustrating a phase offset that 
may be attributable to, for example, Doppler error. 

FIG . 3 is a schematic diagram illustration of one em- 
bodiment of a signal burst or transmission burst, such 
35 as for GSM. 

FIG. 4 is a flowchart illustrating one embodiment of 
a method of the invention. 

FIG. 5 is a table illustrating calculations that may be 
performed by an embodiment of a method of the inven- 
40 tion. 

FIG . 6 is a block diagram illustrating an embodiment 
of a Viterbi traceback reconstructed signal sample index 
in accordance with the invention that may, for example, 
be employed by an embodiment of the invention. 

^5 FIG. 7 is a table illustrating automatic frequency cor- 
rection (AFC) calculations that may be performed by an 
embodiment of a method of the invention. 

FIG. 8 is a table illustrating. a sequence of Viterbi 
operations in parallel with pipelined digital signal proc- 

50 essor (DSP) operations that may be performed by an 
embodiment of a method of the invention. 

Detailed Description 

55 FIG. 3 illustrates a signal burst or transmission 
burst, such as may be employed in a wireless commu- 
nications system, although the invention is not limited in 
scope to a signal burst having this particular format or 
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structure. As illustrated in FIG. 3, the signal burst or 
transmission burst illustrated comprises a predeter- 
mined number of binary digital signals or bits. In this par- 
ticular embodiment, each frame includes, in succession, 
a series of successive predetermined starting bits, a 
predetermined number of successive binary digital sig- 
nals to be transmitted, a series of successive predeter- 
mined training bits, a second predetermined number of 
successive binary digital signals to be transmitted, and 
a series of successive predetermined ending bits. For 
GSM, for example, there are three starting and three 
ending bits, 58 bits in both portions of the signal burst 
comprising binary digital signals to be transmitted, and 
26 training bits, referred to in this context as the "midam- 
ble," for a total of 148 bits per frame or signal burst. Of 
course, the invention is not restricted in scope to GSM. 

As is well-known, GSM uses a form of signal mod- 
ulation in the baseband known as Gaussian Minimum 
Phase Shift Keying (GMSK). GMSK is described in 
more detail in Digital Phase Modulation , by J. B. Ander- 
son, T. Aulin and C. E. Sundburg, 1986, available from 
Plenum, although, of course, the invention is not restrict- 
ed in scope to GMSK or even to MSK. For example, in 
IS54, an alternate baseband modulation scheme is em- 
ployed. In such baseband modulation schemes, such 
as MSK or GMSK, the bit or binary digital signal stream 
to be transmitted, such as a signal burst, is differentially 
encoded, e.g., baseband modulated to produce a posi- 
tive or negative phase shift representing one or more 
binary digital signals in the signal burst being transmit- 
ted. As previously described, this phase shift modulated 
baseband signal may then be applied to a radio frequen- 
cy (RF) carrier for transmission via a wireless medium. 
Therefore, at the receiving end of the communications 
system, after downconversion and signal sampling, the 
binary digital signals in the signal burst being transmit- 
ted may be obtained by (1) a process, referred to as 
"derotation" in this context, applied to each signal sam- 
ple in the signal burst and (2) then passing the derotated 
signal sample through a minimum least squares error 
(MLSE) equalizer In this context, the term "differential 
encoded digital symbol" refers to a complex signal or 
signal sample at the receiving end of the communica- 
tions system. The binary digital signals to be sent are 
transmitted as an analog signal via the modulation 
scheme employed in the baseband. The analog signal 
is then sampled at the receiving end in the baseband to 
provide the complex signal or signal sample. Depending 
on the modulation scheme employed, a symbol to be 
transmitted may comprise a predetermined set of one 
or more binary digital signals or bits. Furthermore, re- 
garding "derotation," for GMSK, for example, a rotation 
of 90° may be applied to each differentially encoded 
symbol or signal sample in the signal burst transmitted 
via the wireless medium, such as by signal multiplication 
in the baseband of each differentially encoded symbol 
or signal sample by 



where k=0, 1, 2, 3 Of course, the invention is not 

5 restricted in scope to a signaling scheme employing a 
particular direction of rotation or derotation. The direc- 
tion will depend, at least in part, on the particular signal 
modulation scheme employed. Likewise, the phase shift 
applied to "derotate" the baseband signal will depend 

to on the modulation scheme employed. Nonetheless, as 
previously indicated, the Doppler effect due to the rela- 
tive motion between the transmitting and receiving sta- 
tion, such as between a mobile station and a base sta- 
tion, may result in phase rotation error in the received 

'5 signal in comparison with the signal transmitted. This is 
illustrated schematically in the Inphase-Quadrature (I- 
Q) plane in FIG. 2. 

As illustrated, a Doppler shift may result in a fre- 
quency error in the received signal that may translate 

20 into a phase offset error in the complex signal obtained 
at the receiving end of the communications system. This 
offset error may also appear in the derotated signal ob- 
tained as well. In FIG. 2, the phase offset error ^denot- 
ed by the phase difference, d9, between Z and Z. Z de- 

25 notes the actual signaj^transmitted including a Doppler 
phase shift error and Z denotes the^ Doppler corrected 
received signal. Likewise, Z' and Z' denote these re- 
spective signals after derotation in the l-Q plane. 

Various approaches to compensating the received 

30 signals for this phase offset error are known, such as 
described in "Two Stage Doppler Phase Corrected 
TCM/DMPSKfor Shadowed Mobile Satellite Channels, 
" by P. J. Mehane, appearing in IEEE Trans, on Com- 
munications , vol. 41, No. 8, August 1993, herein incor- 

3$ porated by reference. However, as previously indicated, 
such approaches are time-consuming and may also 
"bottleneck" significant signal processing capability in 
an environment having limited resources. A method of 
compensating for Doppler error in a wireless communi- 

40 cations system in accordance with the invention to re- 
duce such bottlenecks involves parallel and pipelined 
signal processing utilizing the computational resources 
available, such as, for example, a digital signal proces- 
sor (DSP). For a method of compensating for Doppler 

4 5 error in a wireless communications system in accord- 
ance with the invention, portions of the signal process- 
ing to be performed may be segmented and "offloaded" 
to another processor or coprocessor. These segmented 
portions of the signal processing may be performed in 

50 advance by the coprocessor in a parallel fashion while 
the digital signal processor is also performing pipelined 
signal processing substantially in tandem. In this con- 
text, the term "pipelining" or "pipelined signal process- 
ing" refers to signal processing performed in a predeter- 

55 mined number of separate segments or stages. In such 
"pipelined signal processing," the processing result of a 
particular stage or segment is employed in the process- 
ing performed by the next stage or segment after the 



particular stage. The details of this approach will be- 
come clear in the discussion that follows. 

This approach may be illustrated at least in part by 
the block diagram in FIG. 1. In this particular embodi- 
ment, a digital signal processor (DSP) 170 has embed- 
ded within it a Viterbi decoder 110. An example of such 
a digital signal processor is the DSP1 61 8 available from 
AT&T Corp., which includes an embedded error correc- 
tion coprocessor (ECCP) operating as Viterbi decoder 
110, described in the preliminary data sheet, dated Feb- 
ruary 1994, available from AT&T Corp,, herein incorpo- 
rated by reference, although the scope of the invention 
is not limited in this respect. As illustrated, DSP 170 ob- 
tains the received signal burst. It will, of course, be ap- 
preciated that some preprocessing has typically been 
performed on the signal burst, such as downconversion, 
analog-to-digital conversion, and u de rotation." Thus, 
each symbol in the received signal burst takes the form 
of a complex digital signal in the Inphase-Quadrature (I- 
Q) plane that itself represents one or more binary digital 
signals being transmitted, as previously described. 

As illustrated, and as is well-known, the signal burst 
may be applied to a Viterbi decoder, such as Viterbi de- 
coder 1 1 0, one signal sample at a time in order to obtain 
the transmitted signal burst based on the received signal 
burst. Viterbi decoding is well-known and described in, 
for example, Digital Communications , by E. Lee and D. 
Messerschmitt, available from Kluwer Academic Pub- 
lishers, 1992, Digital Communications by Satellite , by 
Bhargava, Haccoun, Matyas, and Nuspl, available from 
John Wiley & Sons, Inc., 1981, and Digital Communica- 
tions by Satellite , by J. J. Spilker, Jr., available from 
Prentice-Hall, Inc., 1977, all of which are herein incor- 
porated by reference. Viterbi decoding is likewise de- 
scribed in "Maximum Likelihood Sequence Detection in 
the Presence of Intersymbol Interference," by G. D. For- 
ney, Jr., and available in IEEE Trans, on Information 
Theory , IT-18(3): 363-378, May, 1972, and "The Viterbi 
Algorithm," IEEE Proceedings , March, 1973, 268-278, 
herein incorporated by reference. As illustrated in FIG. 
1 , however, prior to being applied to Viterbi decoder 1 1 0, 
each signal sample is adjusted by a phase offset pro- 
vided by Doppler phase 150. Likewise, as illustrated in 
FIG. 1, this phase offset is obtained by Doppler phase 
150 based on a prior predetermined-sized grouping of 
signal samples of the received signal burst that has 
been delayed by time delay 160 and compared with a 
signal estimate of the symbols transmitted for that prior 
grouping provided by Viterbi traceback reconstructed 
signal sample index 140. Delay 160 is introduced to en- 
sure that the appropriate grouping of received signal 
samples is compared with the appropriate signal sample 
estimates. Channel estimate 120 provides signals to 
signal sample index 140 so that signal reconstruction 
may be performed. Signal reconstruction may be per- 
formed in accordance with any method known. Signal 
sample index 140 obtains- signals from Viterbi decoder 
110 based on a process performed by the decoder des- 



ignated in FIG. 1 as a "short traceback" or "short decode. 
0 In this particular embodiment, a short traceback refers 
to a traceback of length one, although short tracebacks 
of length greater than one may also be employed, such 

5 as described, for example, in aforementioned patent-ap- 
plication serial no. 08/152531 , entitled "Variable Length 
Tracebacks." Likewise, as illustrated in FIG. 1, Viterbi 
decoder also provides signals based on a process re- 
ferred to as a "long traceback" or "long decode." Thus, 

10 Viterbi decoder 1 1 0 has the capability to decode the re- 
ceived signal burst in order to determine the binary dig- 
ital signals that have been transmitted. However, a Vi- 
terbi decoder for this embodiment of a method of com- 
pensating for Doppler error in a wireless communica- 

*5 tions system has the ability to perform at least two proc- 
esses, one referred to as a "long traceback" and another 
referred to as a "short traceback." 

One aspect of a method of compensating for Dop- 
pler error in a wireless communications system is a par- 

20 allel and pipelined processing approach between Viterbi 
decoder 110 and DSP 170 to reduce bottlenecks that 
DSP 170 might typically encounter during signal 
processing. More specifically, and as described in more 
detail hereinafter, in processing a signal burst, such as 

25 illustrated in FIG. 3, for example, for a predetermined- 
sized grouping of the received differentially encoded 
symbols, referred to in this context as an intermediate 
grouping, in this particular embodiment each signal 
sample is provided to Viterbi decoder 1 1 0 in order to per- 

30 form a Viterbi update and short traceback. Furthermore, 
in parallel with that process, for the immediately suc- 
ceeding grouping of complex signal samples in the sig- 
nal burst, referred to in this context as the second group- 
ing, digital signal processor 170 forms by pipelined 

35 processing an estimate of the Doppler error in accord- 
ance with the short traceback previously performed for 
each signal sample in the grouping preceding the group- 
ing currently being applied to the Viterbi decoder, re- 
ferred to in this context as the first grouping. Likewise, 

40 the digital signal processor will adjust or correct each 
complex signal sample in the succeeding or second 
grouping in accordance with the estimated Doppler er- 
ror. This is illustrated in FIG. 1, for example, in which a 
phase offset is applied to the received signal corre- 

4 5 sponding to the complex signals or signal samples pre- 
viously described before the signal samples are applied 
to the Viterbi decoder. Thus, for this particular embodi- 
ment, at a given time, an intermediate grouping may be 
processed by the Viterbi decoder, a first grouping just 

50 processed by the Viterbi decoder may be processed by 
the DSP to estimate Doppler error and a second group- 
ing about to be applied to the Viterbi decoder may be 
adjusted based on the Doppler error estimated from 
processing the first grouping. In this fashion, in this par- 

55 ticular embodiment, the digital signal processor typically 
does not experience a bottleneck because processing 
by the digital signal processor continually processes in 
a pipelined fashion complex signal samples based on 
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complex signal samples previously processed in paral- 
lel by the Viterbi decoder, while the Viterbi decoder in 
parallel continually processes additional complex signal 
samples that have been previously adjusted in phase 
by the DSP. 

In a method of compensating for Doppler error in a 
wireless communications system in accordance with the 
invention, processing is applied to the received signal 
burst in advance of the pipelining previously described. 
For example, as previously indicated, the received burst 
signal is downconverted to provide a baseband signal 
and the baseband signal typically is converted from an 
analog signal to a complex signal in a quantized binary 
form, although the scope of the invention is not limited 
in this respect. This quantized binary digital signal rep- 
resents a complex signal in the inphase-quadrature (I- 
Q) plane corresponding to one symbol of a plurality of 
differentially encoded symbols for a signal burst. Thus, 
the received signal burst comprises a plurality of com- 
plex signals or signal samples transmitted via a wireless 
communications system and may be provided to a proc- 
essor, such as digital signal processor 1 70, in the signal 
form just described. As previously indicated, the trans- 
mitted signal burst has a substantially predetermined 
structure in which a subset of the binary digital signals 
being transmitted are known at the receiving end of the 
wireless communications system. As is well-known, this 
signal information may be employed to obtain an esti- 
mate of the communications channel and this channel 
estimate may then be used in further signal processing. 
In order to obtain this channel estimate, the received 
complex signal samples for the encoded symbols in the 
signal burst may first be "derotated" and, likewise, an 
automatic frequency correction (AFC), may be applied 
to compensate, for example, for frequency offset error 
attributable to the oscillator employed to downconvert 
the transmitted signal to a baseband signal, although 
the scope of the invention is not limited in this respect. 

Channel estimation is well-known, such as de- 
scribed in "Design and Performance of Synchronization 
Techniques and Viterbi Adaptive Equalizers for Narrow- 
band TDMA Mobile Radio," by G. D'Aria and V. Zingarel- 
li, published in Nordic Seminar on Digital Land Mobile 
Radio Communication , 3rd Proceeding, Sept. 12-1 5th, 
1988, Copenhagen, herein incorporated by reference, 
and a variety of signal processing techniques, such as 
digital signal processing, may be employed. Typically, 
the training bits in the signal burst being transmitted, 
may be employed to obtain an estimate of the commu- 
nications channel. Likewise, "channel windowing" may 
be employed to determine the "maximum energy" por- 
tion of the signal and normalization and signal scaling 
may, likewise, be employed, although the invention is 
not limited in scope in this respect. Once channel esti- 
mation has been performed, as illustrated in FIG. 1, this 
channel estimate may be employed to perform signal 
reconstruction. Any known technique for signal recon- 
struction may be employed. In signal reconstruction, the 



channel estimate is employed to obtain an estimate of 
- the complex signal sample transmitted for each encod- 
ed symbol capable of being transmitted. This may be 
obtained, for example, from the dot product of the chan- 
s nel estimate with vectors corresponding to the Viterbi 
states for the particular modulation scheme employed. 

In this particular embodiment of a method of com- 
pensating for Doppler error in a wireless communica- 
tions system in accordance with the invention, signal re- 
10 construction is based at least in part on the assumption 
that the complex signal samples obtained at the receiv- 
ing end of the wireless communications system corre- 
spond to locally maximum likelihood symbols also de- 
termined at the receiving end of the communications 
is system. Therefore, reconstructed signals may then be 
employed in signal sample index 140, in conjunction 
with signals obtained from a short traceback performed 
by Viterbi decoder 110, to obtain an estimated signal 
provided to Doppler phase 150 for the transmitted sym- 
20 bol. For a predetermined-sized grouping of signal sam- 
ples in the signal burst, this estimated signal or signal 
sample for each symbol may be compared with the re- 
ceived signal sample for the symbol in order to obtain 
an estimate of the phase offset, as described in more 
2S detail hereinafter. 

As illustrated in FIG. 3, for this particular embodi- 
ment the signal burst may be divided or segregated into 
starting bits, ending bits, training bits and the binary dig- 
ital signals to be transmitted. As illustrated in FIG. 1 , the 
30 received signal samples corresponding to symbols are 
adjusted prior to being applied to Viterbi decoder 110. 
Therefore, initialization of the signal adjustments to be 
applied is desirable. In this particular embodiment of a 
method of compensating for Doppler error in a wireless 
55 communications system in accordance with the inven- 
tion, it is desirable to use the training bits in order to in- 
itialize the Doppler portion of the processing for a variety 
of reasons. For example, the binary digital signals trans- 
mitted are known for the training bit portion of the signal 
to burst. Thus, initialization using this portion of the signal 
burst provides a relatively high degree of reliability in 
terms of estimating Doppler error. 

Although the scope of the invention is not limited in 
this respect, typically in wireless communications sys- 
45 terns employing Viterbi decoding, the signal burst 
processing is divided into a forward portion and a back- 
ward portion, referred to in this context as "forward 
equalization" and "backward equalization." For this par- 
ticular embodiment, signal processing is first applied to 
50 the "forward" portion, ranging from bit positions 77 to 
147 in the signal burst illustrated in FIG. 3, for example. 
In FIG. 3, as illustrated, the initial bit position is num- 
bered 0 and the remaining bit positions increase con- 
secutively. One skilled in the art will nonetheless appre- 
55 ciate extension of the approach to the "backward" por- 
tion of the signal burst, such as to bits positions 65 to 0. 

FIG. 4 is a flowchart depicting the pipelined and par- 
allel processing approach employed in this particular 
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embodiment of a method of compensating for Doppler 
error in a wireless communications system in accord- 
ance with the invention. As illustrated in FIG. 4, it is as- 
sumed that "forward" equalization is performed first. 
Nonetheless, the invention is not restricted in scope in 
this respect and, alternatively, backward equalization 
may be performed first. Next, as illustrated in FIG. 4 and 
described below, AI ACC and AQ ACC are initialized using 
the training bits. As previously Indicated, this is a useful 
aspect of a method of compensating for Doppler error 
in accordance with the invention in that initialization may 
be performed using transmitted binary digital signals 
that are known at the receiving end of the communica- 
tions system. 

In this particular embodiment, using the channel es- 
timate previously obtained, the binary digital signals for 
positions 77-86 in the received signals burst are em- 
ployed to obtain estimated signals or signal samples for 
the transmitted symbols and these estimated signal 
samples are compared with the signal samples actually 
received for these symbols to obtain the difference be- 
tween respective inphase components and quadrature 
components of these compared signal samples. It is as- 
sumed in this context that the difference is largely attrib- 
utable to Doppler error. Thus, the channel estimate is 
applied to known bit positions 77-82, 78-83, 79-84, 
80-85, and 81-86 in the training bit portion of the signal 
burst. This provides 5 signal sample estimates corre- 
sponding to 5 successively transmitted symbols. These 
5 signal sample estimates may now be compared with 
the complex signal samples received. The differences 
between the inphase components of the estimated and 
received complex signals are accumulated to obtain a 
AI ACC estimate for this particular grouping of 5 signal 
samples. Likewise, a AQ ACC estimate for this particular 
grouping is also obtained by accumulating the differenc- 
es between the quadrature components, as illustrated 
in FIG. 4. This initialization approach is summarized by 
the first block in the flow chart shown in FIG. 4, In that 
block, Zj represents the received signal samples with 
Doppler phase offset, whereas Shj* represents the com- 
plex conjugate of the reconstructed signal sample 
based upon applying the channel estimate to the train- 
ing bits, as just described. It will be appreciated that al- 
though in this particular embodiment a predetermined- 
sized grouping comprises 5 successive complex signal 
samples corresponding to successively transmitted 
symbols, the invention is not restricted in scope in this 
respect. 

After bit position 86, the initialization portion has 
been completed and as illustrated in FIG. 4, a signal 
processing loop is employed in which Viterbi decoder 
1 1 0 and DSP 1 70 in this particular embodiment perform 
signal processing in parallel so that a pipelined ap- 
proach maybe employed for the signal processing per- 
formed by the DSP, as previously described. In particu- 
lar in this embodiment, as illustrated in FIG. 4, DSP 170 
employs Al ACC and AQ Acc obtained based on process- 



ing of the immediately preceding grouping of complex 
signals to perform Doppler calculations and, likewise, 
using these calculations, performs Doppler signal cor- 
rections of the complex signals in the succeeding group- 
s ing in the signal burst i.e., bit positions 92-96 in this par- 
ticular example for this particular embodiment, before 
those signals are provided to the Viterbi decoder. At the 
same time, the Viterbi decoder is performing Viterbi up- 
dating and short traceback decoding for the grouping 

10 beginning with bit position 87 for this particular embod- 
iment. This Viterbi updating and short traceback decod- 
ing is performed a symbol or signal sample at a time and 
after processing upon each signal sample, the results 
obtained by the Viterbi decoder are employed by DSP 

75 1 70 to compute Al ACC and AQ ACC for the grouping being 
processed by Viterbi decoder 110 in this particular em- 
bodiment. Likewise, in this particular embodiment, the 
Viterbi decoder performs a Viterbi update by performing 
a long traceback. The Viterbi decoder then performs a 

20 short traceback, of length one in this embodiment, as 
illustrated in FIG. 4. Likewise, once processing of this 
grouping by the Viterbi decoder is complete, then the 
DSP may be employed to perform the Doppler calcula- 
tions and signal corrections for another grouping of sig- 

25 nal samples based on the signal processing by the Vi- 
terbi decoder for this grouping while the Viterbi decoder 
is again performing Viterbi updates and short traceback 
decoding for the grouping of signals immediately after 
this grouping. Thus, it will now be appreciated by one 

30 skilled in the art that by employing this parallel process- 
ing structure, the Viterbi decoder continues to perform 
Viterbi updates and short traceback decodes while in 
parallel the DSP operates in a pipeline fashion to per- 
form Doppler calculations and signal corrections for a 

35 grouping of signals based on the recently completed Vi- 
terbi processing of a prior grouping. 

As illustrated in FIG. 4, eventually forward equali- 
zation is complete, in which case Doppler calculations 
and signal corrections are performed for the final group- 

40 ing in the particular signal burst. Likewise, once this is 
complete, the process may be repeated for backward 
equalization. Ultimately, once the backward equaliza- 
tion is performed, the accumulated Doppler error from 
the backward equalization and the forward equalization 

^5 may be combined to perform automatic frequency cor- 
rection (AFC) for the next signal burst to be processed, 
as explained in more detail hereinafter. This automatic 
frequency correction update is performed because the 
Doppler signal corrections may be employed to assist 

50 in tracking phase offset not easily corrected by employ- 
ing other automatic frequency correction techniques. 

Various approaches to the Doppler calculations and 
signal corrections are possible and the invention is not 
restricted in scope to any particular approach. Likewise, 

55 the invention is not restricted in scope to a particular al- 
location of Doppler calculations or signal corrections to 
successive parallel Viteioi operations while processing 
signal samples in a particular grouping. Nonetheless, in 
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this particular embodiment, 5 successive operations of 
Viterbi updating and short traceback decoding are em- 
ployed in each grouping. On each parallel Viterbi oper- 
ation, the DSP performs a different portion of the Dop- 
pler calculations and signal corrections. Likewise, in this 
embodiment, on each parallel operation, once the Viter- 
bi update and short traceback decoding are performed, 
the DSP recomputes the accumulated Al and AG, i.e., 
AI AC q and AQ ACC , for the grouping being processed by 
the Viterbi decoder based on the short traceback decod- 
ing just performed by the decoder. Signal processing by 
parallel Viterbi operations is summarized by the table 
provided in FIG. 5. For example, on the first Viterbi par- 
allel operation, the DSP takes AI ACC and AQ ACC from 
the preceding grouping and obtains an estimate of the 
change in the phase offset from the inverse arctangent. 
Likewise, on the second parallel operation, the DSP ob- 
tains the average phase offset. This average is obtained 
from an accumulation of the phase offsets obtained for 
each grouping in the signal burst processed so far ad- 
justed by the number of groupings processed. Likewise, 
on the third parallel operation, the accumulated phase 
offset is estimated based on the prior accumulated 
phase offset adjusted by a "weighting" function of the 
average differential phase offset. It will, of course, be 
appreciated that this weighting function may also be 
modified adaptively in alternative embodiments or omit- 
ted entirely. Likewise, on the fourth parallel operation, 
the sine and cosine of the estimated phase offset for this 
particular grouping is obtained. Finally, on the fifth par- 
allel operation, Doppler signal corrections are made to 
the next grouping about to be processed by the Viterbi 
decoder such as illustrated in FIG. 1. It will now be ap- 
preciated that this particular approach to pipelining the 
DSP with the Viterbi decoder processing in parallel in- 
troduces a slight lag in the Doppler signal corrections 
performed. It will be appreciated, however, that between 
successive groupings within a signal burst the impact of 
such a lag should not be significant. Furthermore, due 
to the successive nature of the parallel Viterbi opera- 
tions, the effect of any slight lag should remain relatively 
fixed or stable over the signal burst. 

After the Viterbi decoder has performed a Viterbi up- 
date and the short traceback for the next signal sample, 
such as indicated in FIG. 4, signals are provided to the 
DSP by the decoder so that Al and AQ may be obtained 
for the particular signal sample and then Al ACC and 
AQ ACC may be processed by the DSP. More particularly, 
once a short traceback is completed by the Viterbi de- 
coder, an estimate of the complex signal transmitted is 
obtained using a Viterbi traceback reconstructed signal 
sample index in accordance with the invention, as de- 
scribed in more detail hereinafter although the scope of 
a method of compensating for Doppler error in accord- 
ance with the invention is not limited in this respect. This 
estimate may then be provided to the DSP and the DSP 
may process the dot product of the complex conjugate 
of this signal sample estimate for the symbol with the 



actual signal sample for the symbol in order to obtain Al 
and AQ, although the scope of the invention is not limited 
in this respect. Likewise, Al and AQ for each of the signal 
samples obtained, as just described, may be accumu- 
s lated so that once the final signal sample in a particular 
grouping is processed, AI ACC and AQ ACC may be ob- 
tained for that grouping. As previously indicated, a sim- 
ilar approach to processing Al ACC and AQ ACC may be 
employed in conjunction with the training bits for the pur- 

10 poses of initialization. 

FIG. 8 is a table illustrating the sequence of parallel 
Viterbi operations that may be performed by an embod- 
iment of a method of compensating for Doppler error in 
a wireless communication system in accordance with 

7$ the invention. The table in FIG. 8 illustrates the parallel 
Viterbi operations with the pipelined DSP processing for 
bit positions 87 to 92. As illustrated, this embodiment 
employs 5 parallel operations or operation cycles, al- 
though the invention is not limited in scope in this re- 

20 spect. Therefore, on the Viterbi operation performed by 
Viterbi decoder 110, in parallel the DSP performs a dif- 
ferent portion of the Doppler calculations and signal cor- 
rections, as previously described and as illustrated in 
FIG. 8. An advantage of this pipelined parallel process- 
es ing approach is that the utilization of the DSP and the 
Viterbi decoder is improved in comparison with alternate 
approaches. More particularly, less idle time occurs for 
the processor. 

Although the invention is not restricted in scope to 

30 performing automatic frequency correction or to this par- 
ticular technique of automatic frequency correction, in 
one embodiment, after a phase offset, 0(n), has been 
accumulated based upon forward equalization attribut- 
able to Doppler error, denoted in this context as 0p(N) 

3S and, likewise, 9 B (N) has been accumulated based upon 
backward equalization, as illustrated in FiG. 4, these 
phase offsets may be employed to perform automatic 
frequency correction (AFC) for later signal bursts to be 
received. In this particular embodiment, as illustrated by 

40 the table shown in FIG. 7, the phase offset [9 F (N) - 0 B 
(N)] may be accumulated over a predetermined number 
of signal bursts and then employed to update the phase 
compensation employed to perform automatic frequen- 
cy correction. 

4S As indicated previously, Viterbi decoding is well- 
known and employed in a variety of technologies, such 
as in wireless communications. As is well-known, one 
aspect of Viterbi decoding relates to performing Viterbi 
updates or trellis decoding using received signal sam- 

so pies at the receiving end of a communications system. 
As received signal samples are processed by the Viterbi 
decoder, a "survival metric" is obtained for each Viterbi 
state in accordance with a Viterbi add-compare-select 
operation. After a predetermined number of signal sam- 

ss pies have been processed, symbol decoding is per- 
formed in a manner referred to as a "long traceback" or 
"long decode" that reflects the most likely set of symbols 
to have been transmitted based upon the accumulated 
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metrics. 

In some situations, however, it may be desirable to 
perform a process referred to as a "short traceback." In 
this situation, it may be desirable to form a local estimate 
of the likely symbol to have been transmitted based up- s 
on recently available signal information. Although this 
estimate may not be as good as the estimate based up- 
on a "long decode" or "long traceback;" nonetheless, in 
some circumstances, it may be necessary or desirable 
to have an early estimate available for signal processing 10 
purposes before a long traceback may be performed. 
One example of this desirability is in the context of per- 
forming Doppler error calculations and signal correc- 
tions. Typically, in these situations, it may be desirable 
or important to obtain the results of a short traceback in 1 $ 
a timely and efficient manner. Typically, performing a 
process referred to as "traceback packing" as part of the 
traceback is a cumbersome process including relatively 
significant processing complexity. In this context, "trace- 
back packing" refers to the process by which binary dig- 20 
ital signals or bits indicated by the Viterbi decoding proc- 
ess to have been transmitted are concatenated for sig- 
nal processing purposes. 

FIG. 6 is a schematic diagram of embodiment 1000 
of a Viterbi traceback reconstructed signal sample index 25 
in accordance with the invention providing a relatively 
efficient technique for performing traceback packing, 
particularly for a short traceback. As illustrated in FIG. 
1 , embodiment 1 000 includes a counter or state counter 
300, a traceback shift register (TBSR) 600, a signal re- 30 
construction table 800, and a comparator 900 coupled 
in a configuration so as to provide the sign bit to the 
traceback shift register from a comparison of binary dig- 
ital signals performed by comparator 900. As illustrated 
in FIG. 6, state counter 300 is coupled so as to provide 35 
digital signals to traceback shift register 600. Further- 
more, traceback shift register 600 is coupled so as to 
provide digital signals to signal reconstruction table 800, 
via address decoder 700 for the embodiment illustrated 
in FIG. 6. 40 

Viterbi traceback reconstructed signal sample index 
1000 is intended for use in conjunction with a Viterbi de- 
coder. As is well-known, typically in Viterbi decoding, up- 
on processing of another received signal sample, a Vi- 
terbi update is performed in which stored signals, re- **5 
ferred to in this context as "metrics, " are updated. Once 
this Viterbi update has been performed, the stored met- 
rics may be provided to comparator 900, as illustrated 
in FIG. 6, in order to perform a short traceback in ac- 
cordance with the invention. so 

As illustrated in FIG. 6, MUXes 91 0 and 920 respec- 
tively couple registers 915 and 925 to the respective in- 
put ports of comparator 900. Although not explicitly il- 
lustrated in FIG. 6, these MUXes may be coupled, for 
example, to Viterbi decoder 110 or a RAM coupled to ss 
decoder 110 so that the "metrics" processed by Viterbi 
decoder 110 are available for comparison by compara- 
tor 900. In accordance with the Viterbi decoding proc- 



ess, for a particular Viterbi state, two metrics are com- 
pared to determine the most extreme value metric for 
that particular Viterbi state. The extreme metric thereby 
obtained for that state may be stored for later compari- 
son with the extreme metrics similarly obtained for the 
other Viterbi states. An embodiment of a Viterbi recon- 
structed signal samples index in accordance with the in- 
vention provides an efficient technique for performing 
this processing, as described in more detail hereinafter. 

As illustrated in FIG. 6, state counter 300, in re- 
sponse to an externally-derived clock enable signal, 
may successively provide digital signals corresponding 
to the Viterbi states capable of being transmitted by the 
particular communications system. For example, in FIG. 
6, fourViterbi states are illustrated using a two-bit coun- 
ter, although the invention is not restricted in scope in 
this respect. As illustrated in FIG. 6, the counter provides 
these digital signals to traceback shift register 600. At 
substantially the same time, MUXes 910 and 920 pro- 
vide two stored metrics to comparator 900 for the state 
indicated by counter 300, the stored metrics being de- 
rived from processing by Viterbi decoder 110. More par- 
ticularly, in accordance with conventional Viterbi 
processing, these signals being compared represent 
the sum of a calculated branch metric and accumulated 
cost for the Viterbi state indicated by counter 300. Like- 
wise, the sign bit obtained by comparator 900 based on 
the comparison of the two metrics is then stored in reg- 
ister 930. Depending on the extreme value of the metric 
obtained for this state, this sign bit may later be trans- 
ferred to traceback shift register 600 as a "short trace- 
back bit," as explained in more detail hereinafter. 

As illustrated, the sign bit from this comparison also 
operates as enable signal EN1 to MUX 935. Signal EN1 
operates, via MUX 935, to transfer the one extreme sig- 
nal of the two signals respectively stored in registers 915 
and 925 to register 940. Thus, register 940 operates as 
a "survival metric" register. Now, the contents of register 
940, via MUX 920, may be compared with the contents 
of register 400, which stores the current value of the ex- 
treme metric for the Viterbi states processed up to this 
point. Comparator 900 then compares the "survival met- 
ric" value for the current Viterbi state as indicated by 
state counter 300, stored in register 940, with the current 
extreme metric value for the previously processed Viter- 
bi states, stored in register 400. Assuming the "survival 
metric" for the current state is more extreme, i.e., greater 
or smaller depending upon the particular embodiment, 
the sign bit from comparator 900, denoted in FIG. 6 as 
EN2, enables the contents of register 940 to be loaded 
into register 400. Likewise, the contents of counter 300, 
indicating the current Viterbi state, is loaded as a digital 
signal into the portion of traceback shift register 600 de- 
noted 610 in FIG. 6. Likewise, the sign bit stored in reg- 
ister 930 is loaded into the portion of traceback shift reg- 
ister 600 denoted as 620 and becomes the "short trace- 
back bit." It will, of course, be appreciated that this proc- 
ess is repeated for each Viterbi state indicated by state 
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counter 300. Thus, if a more extreme metric is obtained 
later, the contents of register 400 is loaded with that ex- 
treme metric signal value and binary digital signals rep- 
resenting the corresponding Viterbi state and short 
traceback bit are loaded in traceback shift register 600. 
It should now be clear to one of ordinary skill in the art 
that once counter 300 has provided signals to traceback 
shift register 600 corresponding to all the Viterbi states 
capable of being transmitted for the communications 
system, register 400 should contain the signal value for 
the most extreme metric of the Viterbi states and : like- 
wise, traceback shift register 600 should contain binary 
digital signals representing the particular Viterbi state 
corresponding to that metric signal value and the asso- 
ciated short traceback bit. As illustrated in FIG. 6, the 
contents of traceback shift register 600 may now be pro- 
vided to a signal reconstruction table 800 with memory 
location addresses that correspond to the concatena- 
tion of the Viterbi states in the form of a binary digital 
signal with a short traceback bit. The contents of the 
memory locations of signal reconstruction table 800 
may contain the reconstructed signal sample values 
corresponding to the associated Viterbi states including 
a short traceback bit. Likewise, alternatively, as illustrat- 
ed in FiG. 6, address decoder 700 may translate or de- 
code the binary digital signals provided by traceback 
shift register 600 into the memory location address cor- 
responding to the particular reconstructed signal sam- 
ple value.- 

One advantage of a traceback shift register in ac- 
cordance with the invention, such as the embodiment 
illustrated in FIG. 6, is that it provides a means to effi- 
ciently address a table of reconstructed signal sample 
values. By comparing "metrics, tt as previously de- 
scribed, the traceback shift register ultimately contains 
a binary digital signal that represents the binary digital 
signal most likely to have been transmitted via the com- 
munications system based at least in part upon recently 
available signal information. Thus, this binary digital sig- 
nal may now be employed to address a memory storing 
estimated signal samples, i.e., reconstructed signal 
samples, for the binary digital signals capable of being 
transmitted by the particular communications system. In 
this particular embodiment, a 2-bit Viterbi state is con- 
catenated with one short traceback bit assuming, for ex- 
ample, a three tap communications channel estimate is 
employed, although the invention is, of course, not re- 
stricted in scope in this respect. 

It will now be appreciated that a Viterbi traceback 
reconstructed signal sample index in accordance with 
the invention offers advantages of speed, power sav- 
ings and efficiency in comparison with traditional Viterbi 
traceback packing approaches. The desired recon- 
structed signal sample value may be quickly and effi- 
ciently addressed using a traceback shift register in ac- 
cordance with the invention, as previously described. It 
will also now be appreciated that, although the embod- 
iment shown in FIG. 6 includes a comparator in which 



the extreme signal value of two binary digital signals, 
such as the smaller signal, is obtained, alternatively, this 
technique may be employed to find the larger of the two 
binary digital signals. 

s While only certain features of the invention have 
been illustrated and described herein, many modifica- 
tions, substitutions, changes or equivalents will now oc- 
cur to those skilled in the art. It is, therefore, to be un- 
derstood that the appended claims are intended to cover 

to all such modifications and changes as fall within the true 
spirit of the invention. 



Claims 

75 - 

1. A method of compensating for Doppler error in a 
wireless communications system employing Viterbi 
decoding, said method characterized by the steps 
of: 

20 for each signal sample in a first predeter- 

mined-sized grouping of received signal samples, 

performing a parallel Viterbi update and short 
traceback decode; for a second predeter- 

25 mined-sized grouping, 

forming by pipeline processing an estimate of 
the Doppler error in accordance with the paral- 
lel short traceback decoding performed for the 
first grouping, and 

30 - adjusting each signal sample in the second 

grouping in accordance with the estimated 
Doppler error. 

2. The method of claim 1 , where in the step of forming 
35 by pipeline processing an estimate of the Doppler 

error is characterized by forming the estimate by 
pipeline processing in parallel with performing a plu- 
rality of parallel Viterbi update and short traceback 
decodes for signal samples in an intermediate 
40 grouping between the first grouping and the second 
grouping. 

3. The method of claim 2, wherein the second group- 
ing is characterized by a grouping immediately after 

45 the intermediate grouping in the signal burst; 

the intermediate grouping comprising a 
grouping immediately after the first grouping in the 
signal burst. 

50 4. The method of claim 2, 

wherein the second grouping is characterized 
by the grouping immediately preceding the in- 
termediate grouping in the signal burst; 
55 the intermediate grouping comprising a group- 

ing immediately preceding the first grouping in 
the signal burst. 
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5. The method of claim 2, 

wherein the step of performing a parallel Vi- 
terbi update and short traceback decode is charac- 
terized by performing a short traceback decode with 
a Viterbi traceback reconstructed signal sample in- 
dex. 

6. A method of compensating for Doppler error in a 
wireless communications system characterized by 
the steps of: 

receiving a signal burst transmitted via the wire- 
less communications system, the transmitted 
signal burst comprising a first plurality of sym- 
bols including a second plurality of substantially 
predetermined symbols; 
estimating the channel for the wireless commu- 
nications system by processing the portion of 
the received signal burst corresponding to the 
second plurality of symbols; 
forming reconstructed signals from the estimat- 
ed channel; 

for each signal sample in a first predetermined- 
sized grouping of received signal samples in 
the signal burst, 

performing a parallel Viterbi update and 
short traceback decode in order to select a 
reconstructed signal sample, and 
comparing the reconstructed signal sam- 
ple with the corresponding derotated re- 
ceived signal sample; and 

for a second predetermined-sized grouping, 

forming by pipeline processing an esti- 
mate of the Doppler error in accordance with 
the signal sample comparisons for the first 
grouping, and 

adjusting each signal sample in the second 
grouping in accordance with the estimated 
Doppler error. 

7. The method of claim 6, wherein the step of forming 
by pipeline processing an estimate of the Doppler 
error is characterized by forming the estimate by 
pipeline processing in parallel with performing a plu- 
rality of parallel Viterbi update and short traceback 
decodes for signal samples in an intermediate 
grouping between the first grouping and the second 
grouping. 

8. The method of claim 7, wherein the second group- 
ing is characterized by a grouping immediately after 
the intermediate grouping in the signal burst; 

the intermediate grouping comprising a 
grouping immediately after the first grouping in the 
signal burst. 



9. The method of claim 7, 

wherein the second grouping is characterized 
by a grouping immediately preceding the inter- 
5 mediate grouping in the signal burst; 

the intermediate grouping comprising a group- 
ing immediately preceding the first grouping in 
the signal burst, 

10 10. The method of claim 6, 

wherein the step of performing a parallel Vi- 
terbi update and short traceback decode is charac- 
terized by performing a short traceback decode with 
a Viterbi traceback reconstructed signal sample in- 

1$ dex. 

1 1 . The method of claim 6 further characterized by, prior 
to the steps of performing, comparing, forming, and 
adjusting for first and second predetermined-sized 

20 groupings, a step of initializing the estimate of the 
Doppler error is performed using the second plural- 
ity of substantially predetermined symbols. 

12. A Viterbi traceback reconstructed signal sample in- 
25 dex including a state counter, a traceback shift reg- 
ister, a signal reconstruction table, and a compara- 
tor, the signal sample index characterized by: 

the comparator coupled in a configuration so as 
30 to provide the sign bit to said traceback shift 

register from a comparison of binary digital sig- 
nals; 

said state counter being coupled so as to pro- 
vide digital signals to said traceback shift reg- 
35 ister; and 

said traceback shift register being coupled so 
as.to provide digital signals to said signal re- 
construction table. 

*o 13. The Viterbi traceback reconstructed signal sample 
index of claim 1 2, characterized in that 

said configuration includes a temporary trace- 
back bit register coupled to said comparator so 
45 as to receive said sign bit; and 

said temporary traceback bit register being 
coupled to said traceback shift register. 

14. The Viterbi traceback reconstructed signal sample 
50 index of claim 1 3, characterized in that 

said configuration is adapted to load the con- 
■ tents of said state counter and said temporary trace- 
back bit register into said traceback shift register in 
response to an enabling signal provided by said 
55 comparator. 

15. The Viterbi traceback reconstructed signal sample 
index of claim 12, in which said signal reconstruc- 
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tion table is characterized by 



a binary digital signal address decoder coupled 
to a memory; 

said memory including binary digital signals 5 
stored in memory locations addressed by said 
address decoder; and 

said stored binary digital signals comprising re- 
constructed signal samples corresponding to 
predetermined states provided as binary digital '0 
signals to said address decoder by said trace- 
back shift register. 
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FIG. 7 
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